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HIGH SPEED CYLINDRICAL ROLLER BEARING ANALYSIS 


I, INTRODUCTION 

CYBEAN Cylindrical EEaring ANalysis) has been created 
to detail radially loaded, aligned and misaligned cylindrical 
roller bearing performance under a variety of operating conditions. 
Emphasis has been placed on detailing the effects of high speed, 
preload and system thermal coupling. Roller tilt, skew, radial, 
circumferential and axial displacement as well as flange contact 
have been considered, -Variable housing and flexible out-of-round 
outer ring geometries, and both steady state and time transient 
temperature calculations have been enabled. The complete range 
of e las tohydro dynamic contact considerations, employing full 
and partial film conditions have been treated in the computation 
of raceway and flange contacts. 

Volume I of this report describes the models and associated 
mathematics used within CYBEAN, The user is referred to the 
material contained therein for formulation assumptions and algo- 
rithm detail. The material present in this. Volume II, is 
structured to guide the user in the practical and correct imple- 
mentation of CYBEAN. Input and output architectures containing 
guidelines for use and' a sample execution are detailed in the 
matter which follows. 

II. INPUT DATA 

CYBEAN requires the preparation of input data which in 
general describes the bearing geometry, properties of the mater- 
ials used, and specifies the imposed operating conditions. 

With these inputs defined, optional solution procedure control 
may be selected. 
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The input architecture has been formulated to impose mini- 
mal initial demands on the user. Data is segregated into 
CATEGORIES which individually address specific characterization 
of the configuration addressed. Any single data set item 
required by the program falls into one of ten distinctly iden- 
tified subsets or categories (Table 13 • Category "ROLLER", 
for example, contains all roller geometry data. Category "CAG" 
details cage information. Items detailing operating parameters 
such as load and speed are entered into input categories "0PER8" 
and "LOAD". 

All data required by the basic program are accepted in 

free NAMELIST format and default values are hardcoded to mini- 
mize the demands on user judgment. Special input data, required 

when the program options are used, is in 80 column card image 

format . 

In each category, the user has the freedom to specify 
all, part, or none of the data. If an item of data is omitted, 
a default value is assiomed. A list of these default values 
is shown in Table 2. Failure to include basic data, e.g., ring 
radius, load, speed, etc. results in a diagnostic abort message. 

Data comprising a category is specified in free format, 
with the restrictions that (1) column one of any card is not used 
and ‘ (2) all pieces of data in any category are separated by 
commas^. No specific sequence of data is required within each 
category. A minimum of two cards is needed to specify a complete 
set of data within a category. 

^On-,.some computers the comma must follow the last significant 
digr.t of an input variable. It is suggested that this restric- 
tion. be observed to avoid the inconvenience caused by compiler 
peculiarities . 
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As an example, consider the category **IRING’*. Three items 
are needed. Using the nomenclature of Figure 1, they are: 

1. RIG - The groove radius of the inner ring 

2. FLGALI - The flange layback angle of a flanged 

inner ring on the left. side. 

3. FLGARI - The flange layback angle of- a flanged 

inner ring on the right side. 

An example of free format data is illustrated- for the 
category "IRING" in Figure 2 , In this case the user wishes 
to describe the geometry of a flanged inner ring. Three cards 
are required. The first card contains a Dollar Sign ($)^in column 
2 followed by the category name "IRING”. The second card is 
used to specify values of input data. Note that free format 
is used throughout and that all pieces of data are separated 
by a comma. The third card contains a Dollar Sign column 

2 and the word END in coliimns 3 through 5 signifying the end 
of data for the category. Column 1 is never used in specifying 
data or category. 

Each data category is used to describe a particular aspect 
of either program use or the bearing configuration. Categories, 
in turn, must be arranged in the sequence noted in Table 1 before 
they can be used to transmit data to the program. 

The following paragraphs will list, in their proper order, 
all categories and the data they contain. In certain cases 
where, at the user's option, categories can be omitted from 
the set, these options are made clear. Likewise, it is also 
clearly indicated when a category must be included in the set, 
regardless of execution options. .If the user wishes to omit a 
category of data, he must still include the two cards: 

$CATAGORY NAME 
$END 

^Different computers may allow or require a different symbol. 
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CATEGORY 1 - SOLUTION CONTROL PA^METERS 
CATEGORY NAME : SOLV 


CATEGORY DESCRIPTION : 

Computer program CYBEAN uses a Newton- Raphs on 
iterative scheme to compute values for ‘the governing 
equilibrium equation set. The user may wish to over- 
ride existing solution control parameters. Those, 
which are permitted as input by the user are con- 
tained within this category. 

DATA ITEMS: DEFAULT 


ITMAX - maximum number of iterations to be used in the 
Newton- Raphson Iteration scheme. 


NPR - debug output print flag. Allows the user to see cal- 
culated results at intermediate steps of solution. NPR 
may be input with the following values : 


NPR =0 No intermediate output printed 

NPR = 1 Divergence messages, intermediate 

equation residues (see Volume I) and 
roller- raceway loads are printed 


NPR =2 All output contained in NPR=1 plus 

the corrections of the variable value.s 
(see Volume I) as calculated in the 
linear equation solver 

NPR =3 All output given for NPR=2 plus the 
solution algorithifi used in this par- 
ticular execution 

NPR = 4 All output given for NPR==3 plus the 
matrix of partial derivatives, 


4 




AL78P023 


DATA ITEMS DEFAULT 

CONVER - convergence criterion used to halt the iteration 
procedure. Solution is said to be obtained when 

N 

( E EQ. )/N < CONVER Cl) 

i=l ^ 

and 


SCCIEQJ^ - |EQj^"b/jEQ 1^>1< 2 (2) 

i=l ± X 

Here, h is the iteration index and EQ^ is the i-th equation 
residue , 
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CATEGORY 2 - PROGRAM LOGIC 
CATEGORY NAME: LOGIC 


CATEGORY DESCRIPTION : 

Within this category, the user is permitted to specif 
values for logic used in a given program execution. User 
provided values dictate the program options. 

All variables in this category are "logical”, and 
have either of two values, .TRUE, or .FALSE, eg., 

MPROP = .TRUE. 

In many cases, additional data will be required from 
the user as a consequence of selecting a specific program 
option. Descriptions of this extra input are found in the 
section "SPECIAL INPUT DATA", starting on page 20. 

DATA ITEMS: 

PLTRNG - Allows the user to obtain a line printer plot of 
the inner and outer ring profiles along their axial effec- 
tive length. 

PLTROL “ Allows the user to obtain a line printer plot of 
the roller profiles along their axial effective length. 

ECHO - Allows the user to echo check the input data. This 
option invokes routines which print the data immediately 
after it has been read. 

V 

Special Program Option Logic : 

The following seven data items permit the user to 
invoke certain special program options. All items require 


DEFAULT 

F 

F 

F 


the user to include additional data (see "SPECIAL INPUT DATA") . 




DATA ITEMS 
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DEFAULT 


COEF - Allows the User to input the influence coefficients 
for the housing or other outer ring support structure. 

MPROP - Allows the user to input material properties for 
the rings, rollers, cage and housing. 

OVREND - Allows the user to input values of roller and 
ring radii at the last 3 positions across the effective 
length. These values will overwrite values computed in 
the program. 

SYMY - A .TRUE, value indicates that the roller and rings 
are symmetric about their respective y-axis. If SYMY is 
.FALSE., the non-symmetric roller and race profiles must 
be read in. 

E VS Lie “ The program uses a slicing technique (see Volume I) 
to compute the roller-raceway loads. By default all slices 
are of equal width. In many cases it may be advantageous 
to specify that these slices are of unequal width. Unequal 
slice widths may be included as part of the input data by 
specifying EVSLIC = .FALSE.. 

FITS - Allows the user to use the clearance change portion 
of the analysis. Default is no fit calculation, and a 
flexurally rigid outer ring is assumed, 

THERM - A .TRUE, value allows the user to use either the 
steady state or time transient temperature calculating 
routines . 
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CATEGORY 3 - -ROLLER- GEOMETRY -DATA- - 
CATEGORY NAME: ROLLER 


CATEGORY DESCRIPTION 

Within this category the user must describe the 
geometry o£ the rolling elements within the bearing 
complement. This category is always included. 

(See Figures 3 and 4, all lengths are in mm.) 


DATA ITEMS 


DEFAULT 

ROLLD 

- Roller maximum diameter 

None 

RTL - 

Roller total length 

None 

RCR - 

Roller Crown Radius 

None 


SPHR - Roller end sphere radius on the right side of roller 333 ^ 
(see Figure 3) . 

SPHL - Roller end sphere radius on the left side of roller 381 mm 

RFL - Roller flat length None 

ELO - Effective length of the roller outer ring contact^ None 

ELI “ Effective length of the roller inner ring contact^ None 

XL, XR - The x-coordinate of the roller end sphere origin None 
(see Figure- 3) 

^The effective contact length refers to the longest possible 
length which can be used to transmit load between the roller 
and raceway. Typically, this is the roller total length less 
corner radii. -If, however, the -raceway undercuts are excep-. 
tionally large so that the track width is less than the roller 
effective length then the track width should be input. 
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DATA ITEMS 
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DEFAULT 


EPLAYO - Roller end- flange end' play for a -bearing having q_ 

a flanged outer ring (-see Figure 4). 

EPLAYI - Roller end- flange end play for a bearing having .q, 

a flanged inner ring (see Figure 4) . 

DIACL ^ - Bearing diametral clearance (see Figure 4). 0. 

KLUE - Roller geometry flag having the following four pos- 
sible values: 

KLUE = 1; The roller active profile is either fully 

flat or crowned with a flat. Symmetry about 
y is assumed. (KLUE=1 is applicable to the 
roller shown in Figure 3.) 

KLUE = 2; The roller active profile is fully crowned 
and symmetric about y. 

KLUE = 3; The roller active profile is symmetric about 
the y-axis but will be read in by the user 
(see ''SPECIAL INPUT DATA") . 

KLUE = 4; The roller active profile is non-symmetric 
about the y-axis and will be read in by 
the user (see "SPECIAL INPUT DATA") . 

NUMROL - Total number of rollers in the complement. None 

NS - Number of roller raceway slices used in the analysis. 5 

r 

Note: If the laser specifies symmetry about the roller 
y-axis then NS is the number of slices per symmetric half. 

If the user specifies no symmetry then NS is the total 
number of slices. 

^Calculations with out-of-round components require DIACL to 
be calculated.. See Option 2, page 22. 
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DEFAULT 


DHIl - Angular location o£ the ,£ixs.t- element in -t-he • corapTe- 
ment. In the nomenclature o£ Figure 5, the angle is measured 
CCW positive £rom the bearing y-axis. This input is in 
degrees . 
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CATEGORY 4: OUTER RING DESCRIPTION 

CATEGORY NAME: ORING 


CATEGORY DESCRIPTION : 

ORING is used to describe the geometry o£ the 
outer ring. Specification for lobing of the outer 
ring is made through the input of special data and 
will be discussed later. ORING is always included 
with input. Ring geometry is defined in Figure 1. 
Lengths and angles are to be specified in millimeters 
and degrees respectively. 


DATA ITEMS 

ROG “ Groove radius of the outer race 

FLGALO - Flange angle of the flange located on the left 
side of the outer ring 

■ FLGARO - Flange angle of the flange located on the right 
side of the outer ring. If either FLGALO or FLGARO is 
left unspecified, then the outer ring is considered to 
be without flanges. 

DM - Pitch diameter 

KRING - Used to define the geometry class of rings. 

KRING=1; Both raceways have a flat profile (i.e. 

RIG=ROG=0) and ring geometry is symmetric 
about the y~axis. 

KRING=2; Both raceways have a fully crowned profile 
(i.e. .RIG?^0 and ROG?^0) and their geometry 
is symmetric about the y-axis. 


DEFAULT 

0. Cflat) 
No flange 

No flange 

None 

1 
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KRING=3; Both raceway profiles will be read, in as.. 

user specified input, symmetry is assumed 
(see "SPECIAL INPUT DATA", Option 8, page 29) 

KRING=4; Both raceway profiles will be read in by user 
specified input, no symmetry is assumed (see 
"SPECIAL INPUT DATA", Option 10, page 31). 
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CATEGORY 5 - INNER RING DESCRIPTION 
CATEGORY NAME: I RING 


CATEGORY DESCRIPTION : 

Category IRING is used to describe the geometry 
o£ the inner ring. IRING is always included with 
input data. Ring geometry is defined in Figure 1. 


DATA ITEMS 

RIG - Groove radius of the inner ring 

FLGALI , FLGARI - Flange angles . See FLGALO and FLGARO 
in category ORING. 


DEFAULT 

0. (flat) 
No flange 
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I CATEGORY 6 - CAGE DESCRIPTION 
[category NAME: CAG 

CATEGORY DESCRIPTION : 

The data items contained within this category 
are used to describe the geometry of the cage. This 
set of data must always be included. All lengths 
are in mm (see Figure 6). 


DATA ITEMS 

I RIDE - Cage type flag 

IRIDE = 1; the cage is inner ring land riding 

IRIDE =-l; the cage is outer ring land riding 

IRIDE = 0; the cage is rolling element riding 

RLDC - Rail- land diametral clearance 

SRW - Single rail width 

RLD - Rail- land diameter 

CPCLR - Cage pocket radial clearance 


DEFAULT 

+1 

None 

None 

None 

None 
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CATEGORY 7 - OPERATING CONDITIONS 
CATEGORY NAME: OPER8 ' 


CATEGORY DESCRIPTION : 

■ Bearing operating conditions and operating tem- 
peratures are given in this category. Outer ring, 
inner ring and flange temperatures are used to eval- 
uate the properties of the specified lubricant at 
these locations. Bulk temperature (BULKT) is used 
to evaluate the properties of the specified lubricant 
contained in the free space of the bearing cavity. 

This information is subsequently used in the calcu- 
lation of the viscous drag 'force acting upon the 
rolling elements. Data for this category must always 
be included. All temperatures are specified in degrees 
Celcius . 

DATA ITEMS (See Figure 7) 

SS - Shaft speed - RPM 

BULKT - Average temperature of lubricant in bearing cavity 

TRE - Rolling element temperature 

THSG - Housing temperature 

TSHFT - Shaft temperature 

TOR - Outer ring temperature 

TIR - Inner ring temperature- 

TFl through TF4 - Flange temperatures as shown in Figure 7 


DEFAULT 

None 
100 . 
100 . 
100 . 
100 . 
100 . 
100 . 
100 . 
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CATEGORY 8 - LUBRICATION DATA 
CATEGORY NAME: LUBE 


CATEGORY DESCRIPTION : 

Within this category the user specifies lubri- 
cant properties and other data which relate directly 
to the lubricant or to the definition of friction 
related processes. 


DATA ITEMS 

NCODE - The user may specify particular lubricant proper- 
ties or simply select a value of 1 through 4 for NCODE. 

The latter selection obtains lubricant properties from 
a precoded table. Specific values of NCODE and associated 
lubricant properties are shown in Table 3. The user may 
input lubricant properties not in Table 3 by specifying 
NCODE=0. 

ZTO, ZTI - Lubricant replenishment layer thickness ^ at the 
outer and inner rings, respectively, (mm) 

ZTFO, ZTFI - Lubricant replenishment layer thickness at 
the outer and inner flanges, respectively.' (mm) 


^At the present time the magnitudes of the inner and outer reple 
nishment layer thicknesses have not been correlated with flow 
rate, particular lubricants or bearing speed. The user is re- 
quired to establish proper values of the replenishment layer 
thickness. The following guidelines are suggested: 

1) To avoid starvation, replenishment layer thickness 
should be 1 or 2 times the EHD'film thickness. 

2) Because of centrifugal force, intuition suggests the 
outer be thicker than the inner replenishment layer. 


DEFAULT 

1 


7.62x10 

2.54x10 

1.27x10 

1.27x10 
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DATA ITEMS 
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DEFAULT 

5. 


'^CAV - Percent of lubricant occupying the. bearing cavity^ 

0. £ XCAV < 100. 

FP-K - Lubricant friction coefficient, used in the Allen [3]^ *07 

traction model. Typical values lie in the range 
0.05_< FRK £ 0,08. 

AKN “ Computer program CYBEAN uses a model developed by 
Loewenthal [2] to compute EHD film thickness in point and 50. 

line contacts. The term AKN, the lubricant film thickness 
coefficient, appears in that equation. Typical values 
are 18. £ AKN £ SO. 

XMUCG - Coulomb friction coefficient used at the cage po 
pocket-rolling element contact. If ZTO=ZTI=0, XMUCG is .0175 

applied. 

XMURC - Dry coefficient of friction at race contacts. If .0175 
ZTO=ZTI=0, XMURC is applied. 

XMUFL - Dry coefficient of friction at the flange contact. .0175 
If ZTFO=ZTFI=0, XMUFL is applied. 

The following data must be included if NCODE was specified 


as zero: 

VISl - Viscosity of lubricant (CENTISTOKES) at 100° F. None 
VIS2 - Viscosity of lubricant (CENTISTOKES) at 210° F. None 
RHO60 ' - Density of lubricant (gm/cm^) at 15°C. None 
G - Thermal coefficient of expansion (1/C°) None 
COND - Thermal conductivity (watts/M-Deg C] None 


^As with replenishment layer thickness the amount of free lubri- 
cant should be correlated with the operating parameters. At 
this time such correlations do not exist. XCAV values of less 
than 5 percent are recommended. 

^Numbers in brackets designate References listed in Section 6. 
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CATEGORY 9 - BEARING APPLIED LOAD 
CATEGORY NAME: LOAD 


CATEGORY DESCRIPTION ; 

All bearing applied loads, either forces or 
misalignments are specified in this category. This 
data set must always be included (see Figure 5}. 


DATA ITEMS DEFAULT 

FY - Radial load in Y direction (Newtons) 0- 
FZ - Radial load in Z direction (Newtons) 0. 
THETAZ - Misalignment about the z-axis (degrees) 0. 
THETAY - Misalignment about the y-axis (degrees) 0. 
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CATEGORY 10 - SURFACE FINISH AND FATIGUE LIFE DATA 
CATEGORY NAME: LIFE 

CATEGORY DESCRIPTION : 

Self-explanatory. This category must always 
be included. 

DATA ITEMS DEFAULT 

RMSROL - The RMS surface roughness of the roller (microns) .2032 
RMSIR - The .RMS surface roughness of the inner ring (microns)' .254 
RMS OR - The RMS surface roughness of the outer ring (microns) .254 
CIR - Life correction factor^ for the inner ring !• 

COR - Life correction factor^ for the outer ring !• 


^The numbers input for CIR and COR are used to account for 
improved materials by multiplying the raceway fatigue lives as 
calculated by Lundberg-Palmgren . methods . Typical life factor 
values for modern steels are in the range of 2 to 3. 

In the ASME Publication Life Adjustment Factors for Ball and 
Roller Bearings , the Material Factor D and the Material Process 
Factor E should be used multiplicatively as inputs for CIR and 
COR. 

Additionally, if the user is accustomed to using a lubricant 
life multiplier he must also multiply the material factor by 
the maximum lubricant life multiplier or the "multiplier used 
for ideal lubrication". The program will derate for lubricant 
life effects. 
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III. SPECIAL INPUT DATA 

The user activates extended program capabilities by invo- 
king up to a maximum of eleven options. Logic used to activate 
an option was detailed in the preceding material, and is summa- 
rized in Table 4. 

The user may employ as many options as necessary in a single 
program execution. The only restrictions are foimd in the input 
data sequence shown in Table 4 and in the use of specific input 
formats (Appendix A). 

All options require the user to specify additional infor- 
mation. This data follows immediately after the basic catego- 
rized data. 
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III.l OPTION 1: USER SPECIFIF.D MATERIAL PROPERTIES 


The user may specify the material properties of any bearing 
component. This is done by setting the logical variable 
MPROP = .TRUE., The appropriate properties are entered according 
to the card format shown in Figure Al. Unspecified variables 
(i.e., those set to zero or left blank) are assigned the following 
values ; 


Modulus of Elasticity 

Poisson's Ratio 

Coefficient of Thermal 
Expans ion 

Density 


204083 N/mm^ 

. 3 

0.1124 X 10"^ per °C 
7.806 gm/cm^ 
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III. 2 OPTION 2: FIT CALCULATIONS 


This option enables the- user to analyze cylindrical roller 
bearings manufactured with out-of- round outer rings. The user 
is given. complete flexability in specifying the geometry by 
which preload is induced. The four most popular methods for 
generation of the latter with noncircular outer rings are shown 
in Figure 8. Input card format is specified in Figure A2 . 

Variables which describe the shape of the noncircular com- 
ponents follow: 

MEAN RADIUS - This variable is defined in the nomenclature 
of Figure 9 and Figure 10 as 

n _ ^MAX ^MIN 

*^MEAN 4 

It is to be noted that the mean radius must be specified 
for the raceway profile, outer ring- outer surface profile 
and the housing profile (Figures 9 and 10) . 


ECCENTRICITY RATIO - This variable is defined in the nomen- 
clature of Figure 9 and Figure 10 as 


^^MAX ®MIN^ 

This variable is used to describe the magnitude of out-of- 
roundness manufacture into the raceway, ring outer surface 
and housing. 
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LOBE ORIENTATION ANGLE - The lobe orientation angle, 0, 
is ther angle measured clockwise positive from the bearing 
y-axis to the. first lobe (Figures 9 and 10]. Input is in 
degrees . 

NUMBER OF LOBES - Although most bearings which are manufac- 
tured out-of-round are made with 2 lobes, the user may 
input any number of lobes greater than or equal to zero. 

NOTE: Bearing diametral clearance under this option must 

be specified as follows: 

DIACL={(MEAN RADIUS OF OUTER RACEWAY) - (RADIUS OF INNER 
RACEWAY SURFACE) - (ROLLER MAXIMUM ’DIAMETER) }x 2 
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III. 3 OPTION 5: USER SPECIFIED INFLUENCE COEFFICIENTS 

In the fit analysis described in Section III. 2, it was 
assumed that the housing was rigid, therefore, all deformation 
was experienced by the outer ring. In some instances the user 
may wish to treat the outer ring support structure (the housing) 
as a deformable body. In this case, he must supply the "influ- 
ence coefficients"^ for the housing. The influence coefficients 
completely describe the deformation of the housing and are typi- 
cally obtained by using a finite element analysis. 

The user is required to input one coefficient per roller. 
Coefficients are evaluated at each roller location. Note that 
if the coefficients are not input, default is a rigid housing. 
The card image' input format is shown in Figure A3. 


^The influence coefficients, are defined as the outward radial 

deformation of the housing at the i-th roller location due to 
a unit load at roller location 1. (mm/N) 
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HI. 4 OPTION 4: USER SPECIFIED SL ICE WIDTHS 

CYBEAN uses a slicing technique (Ref. 6 § Volume I of this 
report) to compute lubricant traction at the outer and inner 
ring roller contacts. Ordinarily, the slice widths are assumed 
equal, however. In some instances. It may be advantageous for 
the user to specify their Individual and varying extent. It 
may, for example, be desired to obtain greater detail at the 

roller extremities. 

This option is invoiced by specifying. SYMY = .TRUE, and 
EVSLIC = .FALSE. Since symmetry is assumed, the slice widths 
need only be input for a symmetric half of the roller. The 
numbering scheme used is one where the first slice is encountered 
at the roller centerline, the last at the roller end (Figure 
The card format shown in Figure A4 is used for input data. 
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1 1 1 . 5 -OPTION 5: USER INPUT OF SYMMETRIC ROLLER GEOMETRY. 

The specification of roller crown radius and flat length 
is sufficient to describe the geometry of the roller. If the 
user feels this to be insufficient, he may input values for 
the roller radii at specific locations along the roller 
effective length. The input sequence is shown in Figure 11. 

Note that the number of radii to be read in is equal to 
the number of slices plus one. Card formats are given in Figure 
A5 . The specific example of logic used is SYMY = .TRUE, and 
KLUE = 3. . 
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1 1 1 . 6 OPTION 6: OVERWRITE CALCULATED. TROLLER') END RADII - 

WITH USER SUPPLIED VALUES 

In cases where the roller effective length extends for 
nearly the complete roller length, it ^^^ould be convenient for 
the user to specify the last 3 roller radii. In doing so, this 
would accoimt for the blend radius in the specification of the ■ 
roller geometry. Again, note that symmetry is assumed, and 
the user need only specify the last 3 end radii on a symmetric 
half of the roller. 

The format of the input cards is shown in Figure A6 , logic 
used is SYMY = .TRUE, and OVREND =. .TRUE. 

This option’ also invokes Option 9. 
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I I I . 7 OPTION 7: USER SPECIFIED. COMPLETELY VARIABLE 

ROLLER GEOMETRY 

The user may specify the complete detail of roller geometry 
as input by using the options: 

SYMY == .FALSE, and 
KLUE = 4 

When employing this option, the user must specify roller 
radii and slice' widths across both the outer and inner ring 
effective lengths. Note that the number of slices input corres- 
ponds to the total number of roller racev'/ay slices (Figure 12) . 

Card formats for data input are shown in Figure A7. 
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III. 8 OPTION 8: USER SPECIFIED SYMMETRIC RING GEOMETRY 

As with the roller geometry, the user is given the option 
o£ supplying the program with the symmetric ring geometry. The 
ring line o£ symmetry lies in the y-z plane of the- -bearing and 
divides the ring effective length into two equal parts. 

The data required when exercising this option is shown 
in Figure 13, input card formats in Figure A8 , logic used is 
SYMY = .TRUE, and KRING = 3. 
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III. 9 OPTION 9: OVERWRITE CALCULATED CRING) END RADII 

-W-T-T-H- II.S-F.-R -S-TIPPTrT-F.n VALUES' 

This option is similar to Option 6 except that the user 
specifies the last 3 ring radii instead of the last 3 roller 
radii. Note that when Option 6 is invoked, this option is also 
invoked. 

The format of the input data cards is .shown in Figure A9, 
logic used is SYMY = .TRUE, and OVREND = .TRUE. 
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III. 10 OPTION 10; USER SPECIFIED^ COMPLETELY VARIABLE 

RINO OHOMETRY 

The user may specify the ring geometry as input by using 
the options 

SYMY = .FALSE. and 
KRING = 4 

When using this option the user must specify all of the 
ring radii across both effective lengths*. Note also, that the 
number of slices that are input corresponds to the total number 
of roller raceway slices (Figure 14). 

Card formats for data input are shown in Figure AlO. 
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III. 11 nPTTOM lit TEMPERATURE CALCULATIOMS 

CYBEAN may be used to compute either the time transient 
or steady state temperature distribution within a system defined 
by the bearing and its environment. Logic used requires 
THERM = -TRUE. 

The temperature portion of CYBEAN is designed to produce 
temperature maps for an axi symmetric mechanical system of any 
geometrical shape. The mechanical system is first approximated 
by an equivalent system which consists of a number of elements 
having simple geometric shape. Each element is then represented 
by a node point characterized by a mass, surface area, and having 
either a hnown or an unknown 'temperature. The environment sur- 
rounding the system is also represented by one or more nodes. 

With the node points properly selected, heat balance equations 
are formulated by the program for the nodes of unknown temper- 
ature, These equations become non-linear when there is radia- 
tion between two or more of the node points considered. 

The success of the approach depends largely on the realistic 
physical subdivision of the system. ' If the subdivision is too 
fine, there- will be a large number of equations to be solved. 

If the subdivision is too crude, the results are likely to be 
inaccurate . 

The present thermal simulation is restricted to the treat- 
ment of axially symmetric physical systems. .Bearing rings for 
example, fall into this catego-ry and can be represented by an 
element of uniform temperature. For a component or module which 
is not axially symmetric, the user must represent it with an 
equivalent axially symmetric element of approximately the same 
surface area and material volume. 

This section is based upon work presented in [4]. 
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With input -data prepared as described in the following 
paragraphs, CYBEAN will solve the heat-balance equations for 
either the steady state or the time transient conditions and 
produce temperature maps for the physical system. 

INPUT DATA FOR TEMPERATURE CALCULATIONS 

Card formats for data input are listed in Figure All. 

Card 1 

Card 1 is a control card and contains input for both steady 
state and transient thermal analyses. It is riot intended how- 
ever, that both analyses be executed with the same run. 

Item 1 : Highest Node Number' (M3 . The temperature nodes 

must be ^numbered consecutively from one (1) to the highest node 
number. The highest node number must not exceed one hundred ClOO) . 

Item 2 : Number of Unknown Temperature Nodes (N) . It is 

required that all nodes with unknown temperatures be assigned 
the lowest node numbers. The nodes which have known temperatures 
are assigned the highest numbers. 

Item 3 : Common Initial Temperature (TEMP)°C: The temper- 

ature solution iteration scheme requires a starting point, i.e., 
guesses of the equilibrium temperatures. Card 2 allows the user 
to input guesses of individual node temperatures, however, when 
a node is not given a specific initial' temperature, the tempera- 
ture specified as Item 3 of Card 1 is assigned. 
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Item 4 : Pvinch Flag CIPUNCH) : If thePunch Flag is no.t 

zero (0*) of blank, the system equilibrium temperatures along 
with the respective node numbers will be punched according to 
the format of Card 2. This option is useful if, for instance, 
the user makes a steady state run with lubrication, and then 
wishes to use the resultant temperature as the initiation point 
for a transient dry friction run in order to assess the conse- 
quence of lubricant flow termination. 

Item 5 ; "Output Flag" (lUB) . If the "Output Flag" is 
not zero the bearing program output and a temperature map will 
be printed after each call to the bearing solution scheme. This 
printout will allow the user to observe the flow of the solution 
and to note the interactive effects of system temperatures and 
bearing heat generation rates. Since the temperature solution 
is not mathematically coupled to the bearing solution the pos- 
sibility exists that the solution may diverge or oscillate.' 

In such a case, study of the intermediate output produced by 
the "Output Flag" option may provide the user with better initial 
temperature guesses that will effect a steady state solution. 

Two levels of bearing output are permitted. If lUB is 1, the 
rolling element output is not required. If lUB is 2, full 
bearing output is obtained. 

Item 6 : "Maximum Number of Calls to the Bearing Program" 

(ITl). ITl is the limit on the number of Thermal-Bearing iter- 
ations, i.e., the external temperature equilibrium calculation. 
The user must input a non- zero integer such as 5 or 10 in order 
for CYBEAN to iterate to an equilibrium condition. If ITl is 
left blank or set to zero (0) or 1, bearing performance will 
be based on the initially guessed temperatures of the system. 
Temperatures printed will be based on the bearing generated heats . 
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It is lanlikely that an acceptable equilibrium condition will 
be achievecT. However, the temperatures which result may provide 
better initial guesses for a subsequent r;xn than those specified 
by the user. 


ITl also serves as a limit on the transient temperature 
solution scheme, by limiting the number of times the bearing 
solution scheme is called. Each call to the bearing scheme will 
input a new set of bearing heats to the transient temperature 
scheme until a steady state condition is approached or until 
the transient solution time -up limit is reached. 


Item 7 : ’’Absolute Accuracy of Temperatures for the Exter- 

nal Thermal Solution" (EPI). In the steady state thermal solu- 
tion scheme, each calculation of system temperatures occurs 
after a call to the bearing scheme which produces bearing gene- 
rated heats. After the system temperatures have been calculated 
for each iteration, using the internal temperature solution 
scheme, each node temperature is checked against the nodal tem- 
perature at the previous iteration. 


If (^N-1) nodes i then equilibrium 

has been achieved and the iteration process stops. 


Item 8 : "Iteration- Limit for the Internal Thermal 

Solution" (IT2) . After each call to the' bearing program, the 
internal temperature iteration scheme is used to determine the 
steady- state equilibrium temperatures based on the calculated 
set of bearing heat generation rates. If IT2 is left blank 
pr set to zero .(0) , the number of internal iterations is limited 
to twenty (20). 
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Item 9 : "Accuracy for Internal Thermal Solution" (EP2) . 

The use- o-£ -EP-2- is explained in Volume I. If EP2 is left blank 
or set to zero (0), a default value of 0.001 is used. 

Item 10 : "Starting Time" (START) is a time at which the 

transient solution begins, T^; usually set to zero (0). 

Item 11 : "Stopping Time" (STOP) is the time in seconds 

at which the transient solution terminates, T^. The transient 
solution will generate a history of the system performance which 
will encompass a total elapsed time of 

(Tf - T^) seconds 

Item 12 : "Calculation Time Step" (STEPIN) . The transient 

internal solution scheme solves the system of equations (see 
Volume I) : 


t 


k+1 



T 


T = STEPIN 

The user may specify STEPIN. If left blank or set to zero 
(0) , CYBEAN calculates an appropriate value for STEPIN using 
the procedure described in [4]. 

Item 15 : "Time Interval Between Printed Temperature Maps" 

(TTIME) seconds. The user must specify the length of time which 
will elapse between each printing of the temperature map. The 
interval will always be at least as large as the "calculation 
times tep" (STEPIN). 

Item 14 : "Time Interval Between Calls of the Bearing 

Program" (BTIME) . BTIME will always have a value larger than 
or equal to (STEPIN) even if the user inadvertently inputs a 
short-'er interval. Computational time savings result if BTIME 
is greater than STEPIN, however, accuracy might be lost. 
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Card 2 

In the steady state analysis this card is used to input 
initial guesses’ of individual nodal temperatures for unknown 
nodes as well as the constant temperatures for known nodes, such 
as ambient air and/or an oil sump. 

In the transient analysis. Card 2 is used to input the 
nodal temperatures of all nodes at time = , i.e., at the ini- 

tiation of the transient solution. 

Card 5 

With this card, node numbers are assigned to the components 
of the bearing. With this information the proper system temper- 
atures are carried- into the respective bearing analysis.. The 
inner race and inner ring node numbers may or may not be the 
same at the user’s discretion. Similarly, the outer race and 
outer ring node numbers may or may not be the same. The flange 
numbering scheme is shown in Figure 1. 

Card 4 

The bearing analysis accounts for frictional heat generated 
at five locations in the bearing, i.e., the inner race, the 
outer race, between the cage rail and ring land, the bulk lubri- 
cant due to drag and at the flanges. The heat generated at the 
cage-rolling element contact is added to the bulk lubricant. 

This card allows the heat generated to be distributed equally 
to two nodes. For instance, the heat generated at the inner 
race-rolling element contact should be distributed half to the 
rolling element and half to the inner race. The heat developed 
between the cage and inner ring land may be distributed half 
to the inner ring and half to the cage if a cage node has been 
defined otherwise, half to the- bulk’ lubricant . 
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Card 5 

This' card specifies the node numbers and the Heat generation 
rate at those nodes. This card is used to specify where heat 
is generated at a constant rate such as at rubbing seals or 
gear contacts. 

Card 6 

This card type is used to input the numerical values of 
the various heat transfer coefficients which appear in the equa- 
tions for heat transfer by conductivi ty , free convection , 
forced convection , radi at ion and fluid flow . Up to ten coeffi- 
cients of each type may be used. Separate values of each type 
of coefficient are assigned an index number via card 6 and in 
describing heat flow, paths (Card 7 below) it is necessary only 
to list the index number by which heat transfers between node 
pairs . 

Indices- 1-10 are reserved for the conduction coefficient 
A, 11-20' for the free convection parameters, 21-30 for forced 
convection ,' 31-40 for emissivi ty and 41-50 for fluid flow (pro- 
duct of specific heat, density and volume flow rate). 

As an example, for heat transfer by conduction with coeffi- 
cient A of 53.7 watts/M°C one could prepare a card 6 with the 
digit 1 punched in column 10 and the value 53.7 punched in the 
field corresponding to card columns 11-20. If a conduction 
coefficient of 46.7 were applicable for certain other nodes in 
the system one could punch an additional card assigning index 
No. 2 to the value A = 46.7 by punching a ”2" in card column 
10 and 46.7 anywhere within card columns 11-20. 
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Rather than Inputting constant forced convection coeffi- 
cients, optionally, these coefficients can be calculated by the 
program in one of three ways. If the calculation option is 
exercised a pair of cards is used in place of a single card 
containing a fixed value of a,- The contents of the pair of 
cards depends upon which of the three optional methods are used. 

Option 1) a is independent of temperature 'but is calcu- 
lated as a function of the Nusselt number which 
in turn is a function of the Reynolds number 
R , the Prandtl number P as follows,, (cf. [5]): 

w 1C 

a ■ 

Nu - 

where is the lubricant conductivity, L is 

a characteristic length (with the 'units of meters) 
.and K, a and b are constants. 

Option- 2) a is a function only of fluid dynamic viscosity 
and viscosity is temperature dependent. 

d 

a -= c n 

Option 3) a is a function of the Nusselt, Reynolds and 
Prandtl numbers and viscosity is temperature 
dependent. 

Appendix B has been included to aid the user in data pre- 
paration and calculation of heat transfer coefficients. 
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Card 7 

This card defines theheat flow paths between pairs of nodes. 
Every node must be connected to at least one other node, i.e., 
two or more independent node systems may not be solved with a 
single program execution. 

The calculation of heat transfer areas is based on lengths, 
and L2 input- using card 7. Additionally, the type of surface 
for which the area is being calcuated is indicated by the sign 
assigned to the heat transfer coefficient index. If the surface 
is cylindrical or circular the index should be positive, if the 
surface is rectangular the index should be input as a negative 
integer. 

In the case of radiation between concentric axially symme- 
tric bodies, is the radius of the larger body. For radiation 
between two parallel flat surfaces or for conduction between 
nodes, is the distance between them. 

Fluid flow heat transfer accounts for the energy which the 
fluid transports across a node boundary. Along a fluid node 
at which convection is talcing place, the temperature varies. 

The nodal temperature which is output is the average of the 
fluid temperature at the output and input boundaries. If the 
emerging temperature of the fluid is of interest, it is neces- 
sary to have a fluid node at the fluid outlet. At this auxili- 
ary node only fluid flow heat transfer occurs and the fluid 
temperature would be constant throughout the node. Thus the 
true fluid outlet temperature . will be obtained. 

Conduction of heat through a bearing is controlled by index 
51. The actual heat transfer coefficient which contains a con- 
ductivity, area and a path length term is calculated in the 
bearing portion of the program. • The term is based upon an aver- 
age outer race and inner race rolling element contact. 
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Card 8 

This card inputs data required to calculate the heat capa- 
city of each. node in the system. This card type is required 
only for a transient analysis. 
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IV. OUTPUT DATA 

CYBEAN output is structured to present an immediate defi- 
nition of the problem addressed by the intended computation 
prior to detailing the actual calculation results. The perti- 
nent thermal characterization of the system is noted first. 

This is followed by optional line printer plots of raceway and 
roller geometries and a complete category by category itemization 
of the user supplied data set. All default values generated 
within data subsets are displayed. 

The computed data representing the analysis performed is 
presented after the termination of this initial display. The 
first item of design interest is the fatigue life of the bearing. 
Depending upon the input data this value may have been computed 
by incorporation of user supplied values for life modification 
factors. The bearing life display is followed by roller raceway 
contact load distributions and inner ring resultant reactive 
force moment and displacement vectors. 

The solution to be computed, in general, stipulates that 
a given cylindrical rolling element bearing is subject to a 
constant set of speed, load and position vectors. The analysis 
is requested to generate a set of reaction vectors such that 
equilibrium of forces is achieved. A satisfactory solution is 
characterized by the simultaneous individual equilibrium of rings, 
rollers and cage . 

Several computation procedures, and thus several opportu- 
nities, exist for economic calculation of the reaction vectors. 

For example, iteration algorithms may be constructed to elimi- 
nate all user interaction or conversely, integrate the user in 
an,<on line interactive computation m.ode. The first alternative 
reflects a desire to solve all unknowns ’ simultaneously . Equi- 
librium is sought by automated Newton-Raphson iteration of the 
complete field equation set. 
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The computation expense of this approach becomes evident 
when it is recognized that the field equation set addressed 
in this context typically requires inversions of matrices con- 
sisting of 93,000 elements. The second extreme opposite, com- 
putation strategy would depend heavily on operator experience and 
extensive contingency software to generate a rapid solution. 

CYBEAN, in this first edition, has been structured between 
these two extremes. It has been written to ’take advantage of 
user experience which relates to the selection of solution 
starting values but avoids the trap of having computation success 
depend on user ingenuity. Specifically, CYBEAN is structured 
so that multiple executions bracketing a particular operating 
parameter are performed with inspection of intermediate results 
available for the user. Instead of solving for the complete 
set of unknowns, the field equation set is partitioned to seek 
initial equilibrium between the inner ring resultant reacting 
vector and the given bearing operating condition vectors. Having 

obtained ring equilibrium, its position is fixed in space and 
•equilibrium is established for the rollers and cage. 

Equilibrium is always attained for the roller set but not 
necessarily, depending on the asymmetry of the imposed operating 
vectors, simultaneous equilibrium with the rings. 

Characteristically, a solution is obtained which quickly 
defines a reactive operating vector which is not far removed 
from the one initially posed. Resubmission with altered oper- 
ating vectors allows quick design bracketing, gains intermediate 
results, illustrates effects of changes in operating vector com- 
ponents and arrives at specific information more economically 
than would be possible with other methods. A simplified example . 
follows . 
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STAGE 1 

INPUT LOAD VECTOR - FY = +1000 N 

FX = 0 
FZ = 0 

OUTPUT REACTIVE VECTOR- FY = -900 N 

FX = 0 
FZ = 0 

STAGE 2 

INPUT LOAD VECTOR FY = 1200 

FX = 0 
FZ = 0 

OUTPUT REACTIVE VECTOR- FY = -1050 

FX = 0 
FZ = 0 

We are now within 5% of the desired solution which yields 
a reactive vector component of FY = 1000. 

Specific details of the computation algorithm are presented 
in Volume I of this report. 

CYBEAN output has been structured to be self-explanatory. 
The lubricant data evaluated at calculated steady state opera- 
tion is presented first and is then followed by geometric and 
fit information. Cage forces and speed information precede the 
detailed maximum contact stress and lubricant film thickness 
values computed for raceway and flange interactions. 
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Finally, angular roller rotations detailing tilt and skew 
are displayed before a temperature distribution for the complete 
nodal network which simulates the complete system. 

To detail the program capabilities noted above, a sample 
problem was executed using CYBEAN. The program options employed 
are : 

PLTRNG' = .TRUE. 

PETROL = .TRUE. 

MPROP = .TRUE. 

THERM = .TRUE. 

IV. 1 SAMPLE PROBLEM DESCRIPTION 

The cylindrical roller bearing used for demonstration is 
specified in Table 5. The bearing operates at approximately 
2.4 million DN (bore diameter in mm x'rpm). 

The system modelled incorporates a flexurally rigid magne- 
sium housing, an integral steel shaft-spur gear, ambient air, 
oil sump and a dummy ball bearing with the cylindrical roller 
bearing. Elements represented by the 30 'nodes used -in the 
steady state temperature analysis are identified in Table 6. 

IV. 2 SAMPLE PROBLEM OUTPUT DESCRIPTION 

The complete sample problem output is presented in 
Appendix D. • 

Output contained on the first page informs the user of 
the program version and references the most current program 
operation manual. User invoked program options for this execu- 
tion are also noted on this first page. 
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Pages 2 thru 7 are an organized listing of the user speci- 
fied -speGi-a-1 input data nee'ded for system temperature calculations. 
This data completely defines the thermal simulation. 

Pages 8 thru 15 show line- printer plots of the roller and 
raceway active profiles. 

Pages 16 and 17 present an organized list of the basic 
categorized input data. 

Page 18 : The bearing fatigue life as well as individual 

^10 lives of the outer and inner rings are presented. 

The bearing life represents the statistical combination of the 
two raceway lives. The raceway lives in turn reflect the com- 
bined effects of the user input material factors and lubricant 
film thickness factors. Life modification for materials other 
than basic steel is considered. 

The film thickness to surface roughness ratio is used in 
the calculation of the lube life reduction factor. Detailed 
information for this calculation is given in [6] . 

Pages 19-20: The roller- raceway contact loads for the 

i-th roller at the outer and inner rings are defined in the 
R coordinate frame, illiistrated in Figure 16. These forces 
include both elastic and lubricant traction effects. Rollers 
are numbered in ascending order beginning with the roller lying 
on the bearing y-axis and proceeding counter clockwise. 

Page 21: The inner ring applied forces, moments and dis- 

placements constitute system loading information specified by 
the user. The calculated inner ring reactive forces and moments 
result from the vector sum- of all roller- inner ring contact loads. 
Both elastic and friction forces are included. 
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Pages 22 and 23 are the flange induced roller loads. Num- 
bering scheme is identical to the one used for roller raceway 
loads. These loads are defined in the R coordinate frame of 
Figure 16. 

Page 24 displays the sliding speed magnitude at the roller 
end- flange contact. 

The lubricant data shown on the same page is self-explanatory. 
Temperatures at which properties are evaluated correspond to 
the calculated steady state operating conditions. Note that 
the bearing specified in Table 6, used for this example, has 
flanges only on the inner ring. Referring to the output table 
of lubricant properties, although properties at four possible 
flange locations are shown, only those pertaining to the case 
at hand are used in the analysis. 

Page 25: Output contained on this page is for the most 

part, self evident. Note that since the fit analysis was not 
executed outer ring deflections are zero. If the fit option 
were specified, output would show the combined effects of out 
of round and elastic deflection, the sign convection assumed 
is positive inward. 

Cage pocket normal force is that force experienced by the 
roller due t.o interaction in the Z direction with the cage 
web. A negative sign indicates that that roller is pushing 
the cage. 

Epicyclic speeds, printed along with the calculated speeds 
for the user’s reference, are those speeds assumed by the bearing 
components in the absence of slip. These are -useful in assessing 
roller skid. 

Page 26: The Hertzian contact stresses represent maximum 

values for the line (roller- raceway) and point (roller end- flange) 
contacts . 
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Page 27: The lubricant film thicknesses represent minimum 

values -for the lln'e (roller- racewayj and point (roller end- flange) 
contacts. Dry contact was ‘assumed at the flange in this sample 
problem, therefore, films are zero at this contact. 

Page 28: Calculated roller skew and tilt is presented 

to the user as measured in two distinctly different reference 
frames. "Absolute” refers to the rotation the roller experiences 
relative to its initial position. "Relative" refers to the 
rotation the roller experiences relative to the inner ring 
position. The following example illustrates the two conventions. 

Consider a bearing whose four rollers are "frozen" in their 
position of zero absolute skew and tilt. Assume rollers to be 
located at 90° intervals, two being on the Y-axis and two on 
the Z-axis. The inner ring is now rotated about the Z-axis. 

With the ring in its final position, roller still have zero 
"absolute" skew, however, the rollers which lie on the axis of 
rotation appear skewed when viewed from the "relative" reference 
frame of the inner ring. 

Page 29 summarizes the bearing heat generation rates. 

This data is self-explanatory. 

Page 30: Final operating temperatures of the bearing, as 

calculated under' the temperature calculating program option, 
are shown for all nodes. These steady state temperatures are 
in degrees Celcius. 
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V. PROGRAM LIMITATION AND SPECIAL CASES 

CYBEAN is a design- tool. As with any tool, successful 
use requires awareness of intended applicability and inherent 
limitations . 

A. LIMITATIONS 

The user must conform to the following geometric and oper- 
ating restrictions: 

1) The bearing complement may contain no more than 
fifty (SO) rollers. 

2) Flanges may be specified on the outer or inner 
ring, but not on both simultaneously. 

3) Given a cylindrical roller bearing operating 
with specified misalignment, and/or no geometric 
symmetry, one ring (either the inner or outer) 
must be flanged. 

4) This edition of CYBEAN does not accept externally 
imposed axial loads. 

5) Extremely light radial loading^, wherein a single 
roller interacts with the inner ring, will cause 
error termination. 

6) Use of the fit option (i.e. FITS = .TRUE.), requires 
the mean radius of the outer ring outer surface to 
be less than the mean radius of the housing. 


1 


A current estimate used to determine the minimal radial load is 


P . 


min 


C/SO. 


Here, C is the basic (AFBMA) dynamic capacity and P . is the 
.minimum radial load. 
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B. SPECIAL APPLICATIONS 

Some special applications of CYBEAN are: 

1) The user may appro ximate the bearing heat generation 
rate by specifying ITMAX=1. Heat generation rate 
computations under this option are based upon a single 
iteration of the initial guess independent variable 
values . 

The initial guess variable values are obtained by 
solving the governing equation set for the equilibrium 
of elastic forces, 

2) Through the use of the symmetric ring and roller geo- 
metry program options, the user can use the cylindri- 
cal roller bearing program to analyze a single row 
spherical roller bearing. When doing so, ring misa- 
lignment must be set to zero. 

3) The user may, when making several steady state temper- 
ature program executions, use the card punch option 
(IPUNCHt^O) to obtain the temperatures in 80 column 
card format. These provide an economic initial guess 
(see: Temperature Calculations, Card 2) for subsequent 
runs . 
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FIGURE 1 : USER INPUT RING DATA. 
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FIGURE 3: USER INPUT ROLLER GEOMETRY. 






OUTER RING 


^2 RAIL LAND 
DIAMETRAL 




Y 

FIGURE 6 : USER INPUT CAGE GEOMETRY. 





CYBEAN - COMPLIANT OUTER RING 



(NO PRELOAD) 


USER SPECIFIES: 

USER SPECIFIES : 

USER SPECIFIES: 

USER SPECIFIES: 

USER SPECIFIES 

% 

^N 

i'n 

% 

\ 

I 

I 

I 

I 

I 

RO, RR, RH 

RO, RR, RH 

RO, RR, RH 

RO, RR, RH 

RO, RR, RH 


El 

EO 

EO, El 

EO, El 


<I>I 

^0 

<!>o> h 

‘J’o’ ^I 


Nt=2 

N„=2 

N =2 

N =3 


1 

0 

o 

Nj = 2 

o 

Nj=3 



NOMENCLATURE ; 

Rj^ “ radius to neutral axis of outer ring 

I - outer ring cross section moment of inertia 

RO, RR, RH - mean radii of ring O.D. , ring I.D. and housing, respectively^ 

EO, El, EH - eccentricity o£ ring O.D., ring I.D. and housing, respectively 

({•q, - lobe orientation angle of ring O.D. , ring I.D. and housing, respectively 

No> Nj^, Npj - number of lobes on ring O.D., ring I.D. and housing, respectively 

NOTE: IN THESE EXAMPLES THE HOUSING IS ASSUMED CIRCULAR IN PROFILE. 

FIGURE 8: USER SPECIFIED INFORMATION REQUIRED FOR FOUR MOST 

POPULAR MODES OF INDUCING ROLLER PRELOAD. 



•This Profile May Be 

1 ) Raceway P ro f i le 

2) Outer Ring Outer 
Surface Profile 

3) Housing Profile 


FIGURE 9: 



VARIABLES USED IN DESCRIPTION 
OF A 2-LOBED PROFILE. 



This Profile May Be 

1) Raceway Profile 

2) Outer Ring Outer 
Surface Profile 

3) Housing Profile 


FIGURE 10: VARIABLES USED IN DESCRIPTION 


OF A 3-LOBED PROFILE. 



Line of 
Symmetry 

h inner ring 


effective length 



ROLLER RADII 


SLICE WIDTHS 


ROLLER RADII 


FIGURE 11: OPTIONAL SYMMETRIC ROLLER GEOMETRY INPUT DATA 



Inner Ring Effective Length 


^ h Inner Ring 

Effective Length 

1 



FIGURE 12: OPTIONAL ROLLER GEOMETRY INPUT DATA. 





Line of 
Symmetry 


Line of 
Symmetry 


a\ 

4i- 



RADII 


*4 

e 


H inner ring _ 
ffective lengtii 



T 

Y 


RADII 


lib) Inner Ring 


USER NOTE: 

THE NUMBER OF RADII EQUALS 
THE NUMBER OF SLICES PLUS ONE. 


FIGURE 13: OPTIONAL SYMMETRIC RING GEOMETRY INPUT DATA. 



12a) Outer Ring 


12b) Inner Ring 


USER'S NOTE: 

THE NUMBER OF RADII INPUT IS EQUAL TO THE NUMBER OF SLICES PLUS ONE. 

FIGURE 14: OPTIONAL USER INPUT RING GEOMETRY. 




















FLANGED INNER RING 



67 


TABLE 1 



ORGANIZATION 

OF INPUT DATA CATEGORIES 


CATEGORY' 

NUMBER 

CATEGORY 

NAME 

CATEGORY DESCRIPTION 

NECESSARY 

1 

SOLV 

Solution Control Parameters 

No 

2 

LOGIC 

Program Control Logic 

No 

3 

ROLLER 

Roller Geometry Data 

Yes 

4 

ORING 

Outer Ring Description 

Yes 

5 

IRING 

Inner Ring Description 

Yes 

6 

CAG 

Cage Description 

Yes 

7 

OPER8 

Operating Conditions 

Yes 

8 

LUBE 

Lubrication Data 

Yes 

9 

LOAD 

Bearing Applied Loads 

Yes 

10 

LIFE 

•Fatigue Life Data 

Yes 
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TABLE 2 


DEFAULT VALUES FOR UNSPECIFIED VARIABLES 


CATEGORY 

NAME 

VARIABLE 

NAME 

TYPE^ 

DEFAULT I 

VALUE 

CATEGORY 

NAME 

VARIABLE 

NAME 

TYPE^ 

DEFAULT 

VALUE 

SOLV 

ITMAX 

I 

15 

IRING 

RLG 

R 

0.0 (flat) 


NPR 

I 

0 1 


FLGALI 

R 

NO FLANGE 


CONVER 

R 

•1 1 


FLGARI 

R 

NO FLANGE 

LOGIC 

COEF 

L 

. FALSE . 1 

CAG 

IRIDE 

I 

+ 1 (Inner Ring Riding) 


MPROP 

L 

.FALSE. 1 


RLDC 

R 

— 


OVREND 

L 

. FALSE . 1 


SRW 

R 

— 


SYMY 

L 

. TRUE . 


RLD 

R 

- - - 


EVSLIC 

L 

. TRUE . 

OPER8 

SS 

R 

- 


FITS 

L 

. FALSE . 


BULKT 

R 

100“C 


PLTRNG 

L 

. FALSE . 


TRE 

R 

100°C 


PLTROL 

L 

.FALSE. 

I 

THSG 

R 

100"C 


ECHO 

L 

•FALSE. 

1 

TSHFT 

R 

100“C 


THERM 

L 

. FALSE . 


TDR 

R 

100°C 

ROLLER 

ROLLD 

R 

— 

1 

DIR 

R 

100°C 


RTL 

R 

— 


TFl 

R 

100°C 


RCR 

R 

— 


TF2 

R 

100°C 


SPHR 

R 

381.mmCl5 in.) 


TF3 

R 

100°C 


'SPHL 

R 

381.mmC15 in.) 


TF4 

R 

100°C 


RFL 

R 



LUBE 

NCODE 

I 

4 (MIL-L-23699) 


ELO 

R 



ZTO 

R 

7.620x10’^ MM 


ELI 

R 



ZTI 

R 

2.540x10"* MM 


XL 

R 



ZTFO 

R 

1.270xl0'\ MM 


XR 

R 

“ ~ ” 


ZTFI 

R 

1.2 70xl0‘‘* MM 


EPLAYO 

R 

0. 


XCAV 

R 

5 % 


EPLAYI 

R 

0. 


FRK 

R 

0.07 


DIACL 

R 

0. 


AKN 

R 

50.0 


KLUE 

I 

1 


XMUCG 

R 

.0175 


NUMROL 

I 

“ ~ ~ 


XMUFL 

R 

.0175 


NS 

I 

5 






PHIl 

R 

. 0. 

LOAD 

ALL VARIABLES 

R 

0.0 

ORING 

ROG 

R 

0.0 (flat) 

LIFE 

RMSROL 

R 

.2032 MICRONS 


FLGALO 

R 

NO FLANGE 


RMSIR 

R 

.254 MICRONS 


FLGARO 

R 

NO FLANGE 


RMS OR 

R 

.254 MICRONS 


DM 

R 



CIR 

R 

1. 


KRING 

I 

L 

1 

COR 

R 

1. 


^TYPE REFERS TO VARIABLE TYPE, i.e., I = INTEGER VARIABLE, EXAMPLE: ITMAX = 20 

R = REAL VARIABLE, EXAMPLE: CONVER = .01 

L = LOGICAL VARIABLE, EXAMPLE: COEF = .TRUE. 



TABLE 3 




PROPERTIES 

OF FOUR 

LUERI CANTS 




LUBRICANT 

NUMBER 

fNCODE') 

LUBRICANT 

TYPE 

KINEMATIC 
VISCOSITY 
37. 78°C 
(100°F) 
CVISl) 

(cs) 

98. 89°C 
(210°F) 
(VIS23 

DENSITY 

@15.56°C 

(60°F) 

gm/cm^ 

CRH060) 

THERMAL 

CONDUCTIVITY 

W/m/°C 

(cond;) 

THERMAL 
COEFF. OF 
EXPANSION 
1/°G 10"' 
fG) 

FILM 

THICKNESS 

COEFF. 

AKN* 

1 

Mineral Oil 

64.0 

8.0 

0.88 

0.116 

6. 336 


2 

MIL-L-7808G 

17.8 

3.2 

0.95 

0.152 

7.092 

18.2 

3 

Polyphenal 

Ether 

25.4 

4.13 

1.20 

0.119 

7.470 

24.9 

4 

MIL-L-23699 

28.0 

5.1 

1.01 

0.152 

7.452 

18.2 


*Not part of NCODE information. AKN is input separately in- LUBE category. 



TABLE 4 


ORGANIZATION OF SPECIAL INPUT DATA 


SEQUENCE 

PROGRAM OPTION 

LOGIC USED TO INVOKE OPTION 

INPUT CARD FORMAT 
SEE FIG. BELOW 
IN APPENDIX I 

1 

User Input Material 

MPROP = .TRUE. 

A1 


Properties 



2 

Perform Fit Calculations 

FITS = .TRUE. 

A2 

3 

User Input Influence 

COEF = .TRUE, and 

A3 


Coefficients 

FITS = .TRUE. , 


4 

User Input of Slice 

SYMY = .TRUE. and 

A4 


Widths 

EVSLIC = .FALSE. 


5 

ROLLER GEOMETRY 



5a 

User Input of Symmetric 

SYMY = .TRUE. and 

A5 


Roller Geometry 

KLUE = 3 


5b 

Overwrite Calculated 

SYMY = .TRUE. and 

A6 


End Radii 

OVREND = .TRUE. 


■5c 

User Input of All 

SYMY = .FALSE. and 

A7 


Roller Geometry 

KLUE = 4 


6 

RING GEOMETRY 



6a 

User Input of Symmetric 

SYMY = .TRUE. and 

A8 


Ring Geometry 

KRING = 3 ■ 


6b 

Overwrite Calculated End 

SYMY = .TRUE. and 

A9 


Radii 

OVREND = .TRUE. 


6c 

User Input of All 

SYMY = .FALSE, and 

AlO 


Ring Geometry 

KRING = 4 


7 

TEMPERATURE CALCULATIONS 

THERM = .TRUE. 

All 





table 5 


ROLLER BEARING SPECIFICATIONS 


Inner Race 


Bore Dia. 

mm 

(in) 

Raceway Dia. 

mm 

(in) 

Flange Dia. 

mm 

(in) 

Width 

mm 

(in) 

Groove Width 

mm 

(in) 

Flange Angle 



Outer Race 



Outer Dia. 

mm 

(in) 

Raceway Dia. 

mm 

(in) 

Width 

mm 

(in) 

Rollers 



Diameter 


mm 

Length - overall 

mm 

- effective 

mm 

- flat 


mm 

Crown Radius 


mm 

End Radius 


mm- 

Number 




118 (4.6457) 
131.66 (5.1834) 
137.47 (5.4122) 


26.92 

(1. 060) 

14.59 

( . 5746) 

.6 

deg. 

164.49 

(6.4760) 

157.08 

(6.1842) 

23.9 

( . 942) 


(in) 

12.65 

( . 4979) 

(in) 

■14.56 

( . 5733) 

(in) 

13.04 

( .5133) 

(in) 

8.40 

( . 3307) 

(in) 

622.3 

(24.5 ) 

(in) 

381.0 

(15. ) 


10 


Cage 

Land Dia. mm (in) 137.95 

Axial Pocket Clearance mm (in) .020 

Tangential Pocket Clearance mm (in) .'221 

Single Rail Width mm (in) 4.6 


(5.4312) 
( .0008) 
( .0087) 
( .18 ) 


Operating Conditions . 

Shaft Speed 

Bearing Radial Load 4450 N 

Oil Inlet Temperature 366. 5K 


20,000 rpm 
(1000 lb) 
(200°F) 


Misalignment of Races 


15 min. 


Lubricant 
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TABLE 6 


NODE 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 


NODE DICTIONARY 

ITEM 

Shaft 

Input Pinion 
Shaft 

Inner Ring 
Ring Flange #3 
Ring Flange #4 
Rollers 
Outer Ring 
Cage 
Housing 
Housing 
Hous ing 
Shaft 
Hous ing 
Shaft 

Dummy Ball Bearing 

Housing 

Housing 

Housing 

Internal Air 

Lubricant Entering Bearing 
Lubricant in Bearing Cavity 
Lubricant Entering Bearing 
Lubricant Entering Dummy Bearing 
Lubricant Entering Dummy Bearing 
Lubricant in Dummy Bearing Cavity 
Lubricant Entering Dummy Bearing 
Lubricant 
Lubricant in Sump 
Air, External to Housing 
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APPENDIX ,A 

SPECIAL DATA INPUT CARD FORMATS 
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FIGURE A2; .FIT .CALCULATION USER INPUT CARD FORMAT (PAGE 1 of 2) 



NOTES; 1) 'CTBEAN' USES AN ITERATIVE TECHNIQUE TO OBTAIN' THE DEFORMED SHAPE OF THE OUTER RING. THE NUMBER OF- FIT ITERATIONS REFERS TO THIS SCHEME 
IF LEFT BLANK A PRESET MAXIMUM OF 10 ITERATIONS WILL BE PERFORMED.. 

2-) SOLUTION ACCURACY REFERS TO THE CONVERGENCE USED IN (1). IF THE MAXIMUM CHANGE IN DEFORMED SHARE IS l£SS THAN THE SOLUTION ACCUHACY, 
THEN THF rilRBENT DEFORMED SHAPE IS TAKEN TO BE A SOLUTION, if LEFT BLANK A VALUE OF (1/ROLLER DIAMETER) X 10'^ IS ASSUMED. 







FIGURE A2; FIT CALCULATION USER INPUT CARD FORMAT (PAGE 2 of 2) 
CARD 3 



















FIG. A3: USER INPUT OUTER RING INFLUENCE COEFFICIENTS CARD FORMAT 
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FIG.'A4: USER INPUT SLICE WIDTH CARD FORMATS 
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SLICE WIDTH ACROSS SYMMETRIC HALF OF OUTER RING EFFECTIVE LENGTH (USE AS MANY CARDS AS NEEDED, mra) 


WIDTH OF SLICE 
NO. 1 

WIDTH OF SLICE 
NO. 2 

ETC. 








CARD 2 


fl 

B 

a 

a 

a 

a 

a 

a 

fl 

B 

fl 

B 

B 

B 

a 

B 

a 

a 

a 

a 

a 

i 

i 

a 


iS 

a 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

R 

1 

1 

1 

1 

1 

1 

1 

1 

1 

R 

1 

1 

R 

1 

1 

R 

1 

R 

1 

e 

R 

1 

R 

1 

1 

1 

R 

R 

R 

1 

1 

R 

1 

R 

1 

R 

1 

1 

R 

1 

R 

R 

1 

R 

1 


32' 


iaj;3-ijsl36ii7;i6'39:4o;ji 


|j2U3Ui;4S;J6.S7| 


Ualioiji 


5253 


Isj-lsslse 


57 


5e>59i 


!62'63'64,'6566l67i 


|Ee^9 


70 71 


72 73174 175 


761771 


I 71 I 


SLICE WIDTHS ACROSS SYMMETRIC HALF OF INNER RING EFFECTIVE LENGTH 




ETC 


(USE AS MANY CARDS AS NEEPp, min') 




WIDTH OF' SLICE 
NO. 2 

























FIG. A6: USER INPUT END -SLICE ROLLER RADII 

Card 1 



•Card 2 













omimi PAGE IS 
QE quality 


CO 

ts3 


FIG. A7: USER INPUT ROLLER GEOMETRY CARD FORMATS 

CARD 1 


^ 4 


*^7j0 


6IJ7 


^j2C:c:'j22'23l24ias|26'27Vaf59t50i 


32l33j34b5}3637 


3fi‘39| 




~r 


sal 


54;ss 


!S6 


kojol l6a'o3!64!$sii36 [67'eeka ro,7l 


■76 ~rl7«l« 8C 


llo . 0 


1-1 II oi .I 0 


WIDTH OF .SLICE 

WIDTH OF SLICE I 

ETC. 





NO. 1 

NO. 2 1 




- 


0\KD 2 


I a 3!-: S|6|7 8|9 .0,11 12 


16 If? IB |'9 


jaoiai jaaiaslaalasIaelaTtzaiagaj'si 


32'33;3~l3S;36!37i3e39Ue!4 |■^^'43i4J..15|J6;97|Ja{49■£o'5l |i2l£3 |S.^£5 j56|&7|68;S!i;eO,6l kz t3>M;e5]c6!67|g6,63[70_7l |72 Mj?.* TSp -7|78l73[aci 


FI 0 < 0 


F 1 0 . O: 


RADII NO. 1 

RADII NO. 2 

ETC. 






CARD 3 


12 3 4 jB |6 


13 I6il7 


.j ■ ^ r - , ^ ^ — — r — lT..^l 

'^I|^ 23|24|2S!26 27;2a|29‘3S]3l [32;34|34|33;joj37|38'3j|40j4i |4;U3U4:43U6;47j4ai49 Cob |52|t3j54,'55|5lij57[s6|39,eOol i&2 s3;S4'Ei66 !67|eelS9pO 71 172 73j74i7S '^■771 


T” 

|7aiT»laci 


r I 0 . 0 


FI 0.0 


■•SLICi: WTUmS across inner ring effective length ( use as many cards as Nri-PF.n. mml 


W I Di ll Ol- .SLICn 
NO. 1 


WIDTH OF .SMCE 
NO. 2 


ETC. 


CARD <1 


2j3j4j3[6l7l89|lO.,llt2|l3 |l4jl3 


^ j aSiasIzTiaa j as 30j3! i32;33j343S }3c I 37 [3a 3s|4Cl4l bsUa U J :4si46-47 146 i49|M!si 521^3 JB 37 Sa|E9je6oi^'&3l64;6S^6!67;c6^9;70;7f|7g'73i74 7B|76 -7|78|7J 


3 


r 1 0 . 0 


FI 0.0 


ROI.U.I! RADII AT Sl.ICIi I.ND.S ACROSS TIIE.OUTTR RING EFFECTIVE LENGTH (USE AS MANY CARDS AS NEEDED, 


mill) 


RADII NO. 1 


RADII NO. 2 

• ETC. 

- 
























FIG. \.8: USER INPUT SYMMETRIC RING GEOMETRY CARD FORMAT 

CARD 1 

HI" i°i’ 1® i® i'°i" |wjl6|l6jl7jia].j|20'ai'|22’2jia4|25l2fi!2.7|28 


|J~Di I3r: J7 l3a'39|<C|4i |J2 43 U4|45!4e'47l4al49 IsoiSI 


SIS1599Si 


$8 

liJi 

kc 

|6|';52! 

Io3!&4 

■2b| 

iSfij 


74 75i?6»'7!7a|'*9!8C 


SYMMETRIC RING RADII AT SLICE ENDS ACROSS OUTER RING EFFECTIVE LENGTH (USE AS MANY CARDS AS NEEDED, mm) 


RADIUS NO. 1 RADIUS NO. 2 


32 

1 

38i36!37 

i , 

36 

39 

40 


J2'43 


•M|4J'iO 

1 I 1 

51 |62[53 



2l63|&4|65j56|67]60te9i7O,7l |t 2;73 74 7jj76l77l7fl|79 




SYMMETRIC RING RADII AT SLICE ENDS ACROSS THE INNER RING EFFECTIVE LENGTH (USE AS MANY CARDS AS NEEDED; ram) 


RADIUS NO. 1 RADIUS NO. 2 

















FIGURE A9: USER INPUT END SLICE RACEWAY RADII 
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FIGURE All CCONTINUED) : INPUT DATA CAftD FORMATS FOR TEMPERATURE CALCULATION 
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FIGURE All (CONTINUED): INPUT DATA CARD FORMATS FOR TEMPERATURE CALCULATION. 





FIGURE All (CONTINUED): INPUT DATA CARD FORMATS FOR TEMPERATURE CALCULATION, 





















CARD 7 - HEAT FLOW PATHS; USE AS MANY' CARDS AS NEEDED, FOLLOWED BY A BLANK CARD 
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FIGURE All (CONTINUED) : 


INPUT DATA CARD FORMATS FOR TEWERATURE CALCULATION. 


































CARD 8 - NODAL HEAT CAPACITIES (USE ONLY FOR TRANSIENT CALCULATIONS) 


ONE CARD PER NODE ■ 



FIGURE All (CONTINUED): INPUT. DATA CARD FORMATS FOR TEMPERATURE CALCULATION. 








APPENDIX B 


HEAT TRANSFER COMPUTATION NOTES 



B.l 


BASIC EQUATIONS* 


B.1.1 Heat Conduction 

The rate of heat flow q . • (W) rhat is conducted from node 

> J 

i to node j may be expressed by, 




'Cl, 3 


Ct. - t.) 
1 3 


t . and t . are the temperatures at i and j , respectively, A- • the 
^ J 2 ^ > J 

area normal to the heat flow, (m ) L. • the distance (m) and 

the thermal conductivity betx^een i and 3 ', (W/m°C) . 


Assuming that the structure between point i and 3 is com- 
posed of different materials, an equivalent heat conductivity may 
be calculated as follows: 



The calculation of the areas will be discussed in Section B.l. 5. 


*This Appendix is based on the material present in Reference 4. 
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B. 1.2 Convection 


The rate of heat flow that is transferred between a solid 
structure and air by free convection may be expressed by 


‘VI ,1 = a •• 

\ J 


A. : 


"i - 


t . 
3 


1.25 


SIGN (t. - t.) 

3 


where 


( 1, if (ti 

SIGN = < ^ 

if Ctj 


t.) > 0 

r - 

t^) < 0 


in which 


a 


ij 



• 10 ^ W/m^ “ (degC)^’^^ for hot surfaces facing 
upward and cold surfaces facing downxirard 

• 10 ^ W/m^ - CdegC)^’^^ for hot surfaces facing 
downward and cold surfaces facing upward 

‘ 10 ^ W/m^, - fdegC)^'^^ for vertical surfaces 


For other special conditions, a-- must be estimated by referring 
to heat transfer literature. 

The rate of heat flow that is transferred between a solid 
structure and a fluid by forced convection may be expressed by 


q . . = a. -A- . ft. - t-1 
^ni,j '■ 1 j-* 

in which a. . is the convective heat transfer coefficient, 
ij 

Now, with a'= > introduce the Nusselt number 

XT OiL 


the Reynolds number 


\ = 

and the Prandtl number 


UL 


P = P 
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where 


L is a characteristic length which is equal to the diameter in 
the case o£-a cylindrical surface and is equal to the plate 
length in case of a flat surface (in) • 

U is a characteristic velocity which is equal to the difference 
between the fluid velocity at some distance from the surface 
and the surface velocity (m/sec) . 

A is the fluid thermal conductivity (W/M°C) 

2 

V is the fluid kinematic viscosity (M /sec) 

^ is the fluid density (kg/m^) 

Cp is the fluid specific heat (J/kg°C) 

For given values of R and P the Nusselt number N and 
® e r u 

thus the heat transfer coefficient may be estimated from one of 
the following expressions: 

Laminar flow along a flat plate: < 2300 

N = 0.323 /R • ^/P 
u e r 

Laminar flow of a liquid in a pipe: 

N^=1.36 Vr^.P^C^) 

where D is the pipe diameter and L the pipe length 

Turbulent flow of a liquid in a pipe: 

N = 0.027 • R*^'® • ^/P 
u • e r 

Gas flow inside and outside a tube: 

0 5 7 

•*N = 0.3 R 

u e 


Liquid flow outside a tube: 


N, = 0.6 R^*^ 
u e 


,0. 31 
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Forced convection from the outer surface of a rotating 
shaft 

N = 0.11 [ 0.5 R ^ • P 1°-^^ 
u e r"^ 

where the Reynolds number R is developed by the shaft 
2 ® 

_ MTT D 

V 

in which m is the angular velocity [rad/sec) 

D is the roll diameter (m) 

The average coefficient of forced convection to the lubri- 
cating oil within a rolling contact bearing may be approximated 
by, 

a= .0.0986 [1 + }^ACP 

using + for outer ring rotation 
- for inner ring rotation 

in which N is the bearing operating speed (rpm) 

D is the diameter of the rolling elements (mm) 

dj^ is the bearing pitch diameter (mm) 

6 is the bearing contact angle; zero for cylindrical 
roller bearings (degrees) 

B.1.3 Fluid Flow 

The rate of heat flow that is transferred from fluid node 
i to fluid node j by fluid flow is 


lx* • “ pV- . C (t . 

^fi,j 13 p- 1 


tp 


is the volume rate of flow from i to j' . It must be observed 
that the continuity of mass requires the following equation to 


be satisfied 

Et. . = 0 

ij 

provided the fluid density is constant. The summation should be 
extended over all nodes i within the fluid which have heat 
exchange with node j by fluid flow. 
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B.1.4 Heat Radiation 


The rate b£ heat flow that is radiated to node j from 
node i is expressed by 


whe re 


^Ri,j “ '^i,j 


T- = t- + 273.16 

v* J 

T. = t. + 273.16 
1 1 


{(t. + 273)"^ - Ct. + 2 73)^} 

^ J 


and the value of the coefficient. 6 . . depends on the geometry and 
the emissivity or the absorptivity of the bodies involved. 

For radiation between large, parallel and adjacent sur- 
faces of equal area, A- . and emissivity, c. 6- - is obtained 
from the equation 

5 . . = £ . . a A. . 

3 .,] ■ 

where a, the Stefan -Boltzmann constant, is 
a = 5.76 • 10"^ W/m^/CdegK)^ 

For radiation between concentric spheres and coaxial 

cylinders of equal emissivity, s- .,6. -is given by the 
equation 

e • . oA . . 

. = i>3 

•^ij A. . 

l+(l-e- ■) 


where o is as above A. . is the area of the enclosed body and 
A*. . is the area of the surrounding body, i.e., A. . < A*. .. 

■ 1,3 1,3 

Expressions for 6. • that are valid for more complicated 
geometries or for different emissivities may be found in the 
heat transfer literature. 
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B.1,’5 Calculation o£ Areas 


In the case of heat transfer in the axial direction ^ 
is given by the equation (Fig- B-3) 

A. . = 2llr • At 

1, j m 


Referring to the temperature calculation input instructions, 
card 7, but recalling L must be input .in mm not m. 




In the case of heat transfer in the radial direction, 
A. . is obtained from the expression 

1 > j 


A. . = 2Jlr 

m 


H; L, = r^; L. = H 




and similarly for the radiation term above 

A*. . = 2Er* H 
i,j m 

o m 


Lz = 2H 

in which H is the length of the cylindrical surface; where heat 
is conducted between i and j , is given by the same equation 
as above CFig> B-4 (a)); where heat is convected between i and j, 
r^^ is the radius of the cylindrical surface (Fig. B-4(b)); where 
heat is radiated between i and j , is the radius of the enclosed 
cylindrical surface and r^^* the radius of the surrounding cylin- 
drical surface (Fig, B-4(c)). ■ 
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Fig. B-3 


Fig. B-4(a) 


B.2 TRANSIENT ANALYSIS 

For the transient analysis all o£ the data pertaining to 
the node to node heat transfer coefficients must be provided by 
the input. Additionally, the volume and the specific heat at 
each node is required. For metal nodes this input is straight- 
forward. However, when fluid flow is being considered there is 
no easy way to approximate the fluid nodal volume in a free 
space such as the bearing cavity. However, through use of CYBEAN 
the user's ability to make appropriate estimates will improve. 


Fig.B -4(b) 


,Fig. B -4(c) 


KM 




APPENDIX C 
CYBEAN FLOWCHART 


101 




102 





COMPUTE DERIVED QUANTITIES 


HEAT GSflERflTlON RAT3S 


i>i > fatigue life 


MATEEUAL MODULUS EFFECT 


RACEWAY-R.E. CONDOCTIVITy 
FLAMGS -R.E. CONDUCTIVITY 


WRITS RESULTS 


VES DOINGN 

'“'^MP.CALC,?^ 


STEADY STATE 


SOLVE TEMPERATURE DISTRIBUTION 
FOR STEADY STATE EQUI LIBRIUM 


SOLVE TEMPERATURE DISTRIBUTION FOR 
TRANSIENT TIME STEP fi WRITE RESULTS 







APPENDIX D 
SAMPLE OUTPUT 
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Tftu'uaTHas' ifji'. — * — 1 . ' V a t. a"H"7 ~ f r' »" 5 k 


CYflEAN/N^SA - VERSION NUMBER 1 {ORIGINAL) 

L'SE THIS VERSION JITH ThE FOLLOWING USER"S MANUALS ' 

- VOLUME I REVISION 0, DATED S» t H >1% S K F REPORT ND . AUJ @> POZ|3|. 

- VOLUME 2 REVISION 0i DATE 06 S K F REPORT AUT^ P02’3 

program options ineffect 


SYiIHETRIC roller AMO RING GEOMETRIES ABOUT Y-AXIS 
EQUAL SLICE WIDTHS 
PLOT RING profile 
PLOT ROLLER PROFILE 

USER INPUT MATERIAL PROPERTIES"' " — 

TEMPERATURE OI3TRIOUTION WILL BE COMPUTED 
ALL DATA PRESENTED IN METRIC UNITS 


units 


ENGLISH 


METRIC 


length 

FILM THICKNESS 

MASS 

FORCE 

COEF. OF THERMAL" EXP. 

DENSITY 

STRESS 

MO'ICHT 

TEMPERATURE 

PRESS. VIS. COEF. 

SPEED 

HEAT GC-J. rates 
replenish, layer thick. 


INCHES 

INCHES 

SLUGS 

POUNDS 

1/OEGREES F 

LBS. /INCHES CUBED 

PS I 

INCH-POUNDS 

DEGREES F 

LB-INCHES SOUAREO 

RPM 

WATTS 

inches 


MILLIMETERS 

MICRONS 

KILOGRAMS 

Nf.ViTONS 

1/DE6R£ES"C' "" ' 
NTS./MH CUBED 
NTS./M.M SQUARED 
MH-NEWTONS 
DEGREES C 
NT-MH SOUARED 
RPM ■ • 

WATTS 
MICRONS 


Samnle Outiout - Page 1 
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Tl TOUS T HTI 'rTTrr: — I c Y 'B E A N 7 W A' S » 


3TCAUY GTATC TCMPCRATUHK CALCULATION. ITEIIATION LIMITS 5 
NODE POINTERS 

O.RACE I. PACE BULK OIL FLNG.l FLNG.2 FLHG.i 
8 A 22 8 8 5 


ABSOLUTE ACCURACY = 2.O0OE5REES 


FLNG.A cage SHAFT 
,6 0 3 


l.RINO ROLL. EL. O.RING HSG, 

•4 -J g. . 


NODES WHERE BEARING HEAT IS GENERATED 

OUTER RACE INNER RACE R.E.DRAG 

T- 8 4-7 4-11 


CAGE-R.C. 

t-4 


CAGE" LAND FLNG.I-RE 

5-9 7-& 


■FLNG.2-RE" FLNG:'3“RE FCa1g;<i-rc' 

7-& 5^-7 6-7 


CONSTANT GENERATE HEATS 

node gen. heat 


NOOC GEN. HEAT 


NODE "GEN. HEAT ' "‘NODE GEN. HEAT NODE 'G‘EN‘. H.EaT'' 


2 300.00 


IS ' 300.00 


16 lAOO.OO ■ 17‘ 300.00 


S^ple Output - _Page 2 
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' * ■ " V » TL A " ! i * ' '"r'mHNumu i' u ' lv ' i^ TCi? "'aKf’'’"'~nTt)TnnKit.i>"'tm:. > i^Tn^AVi/MAiiA 


heat TRA-JSFFS COEFFICirtJTS 

TYPE IN3EX 


COEFFICIENTS 


COHOUCTION 
COMOUCTION 
CONOUCTIOY 
FORCED CONVECTION 
FORCED CONVECTION 
FORCED CONVECTION 
FLUri) flo-j 
FLUID FLO« 


1 

A6.7000 

2 

53.7000 

5 

SO. 0000 

21 

3000.00 

22 

'109.100 

23 

37.0(100 

R1 

50.5000 

<l2 

119.000 


C ^ 

^ s> 

*< Sr 


- Page 3 
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** * C Y B E A « / M A S V"* 'TECHNOUOGr DIVISION ~SKF INDUSTR IES“ INC C“ Y" B "E" A “IT /~N"A'-S~S 


DESCRIPTION OF THE GEOMETRY AND INDICATION OF THE TYPES AND PATHS OF HEAT TRANSFER 


ALL lengths are ‘IN 

MILLIMETERS, A 

NEGATIVE 

SIGN OF 

THE INDEX MEANS 

NO ROTATIONAL 

SYMMETRY 

TYPE OF HEAT TR . 

INOCX 


NODE 


NODE 

1ST LENGTH 

2ND LENGTH 

3R0 LENGTH 

CONOUCTIOH 

i 

BETWEEN 

1 

AND 

2 

20.0000 

59.0000 

"BICOOO"* 

CONDUCTION 

1 

BETWEEN 

1 

and 

3 

R«0000 

55.0000 

33.5000 

CONDUCTION 

1 

BETWEEN 

3 

and 

A 

26.9200 

59.0000 

7.5000 

CONDUCTION 

1 

BETWEEN 


AND 

5 

6.8300 

62.AOOQ 

“lO.AOUU 

CONDUCTION 

1 

BETWEEN 

4 

AND 

0 

6.8300 

62.A000 

lO.AOOO 

COMQUCTIO *i 

1 

aCTWEEN 

3 

AND 

13 

ft.OOOO 

55.0000 

46*5000 

CONDUCTION 

1 

BETWEEN 

13 

AND 

15 

a*0Q00 

55«0000 

A3.C0D0 

CONDUCTION 

1 

BETWEEN 

15 

AND 

16 

1.0000 

0.1091 

16.00UU 

CONDUCTION 

2 

between 

10 

AND 

lA 

5. 0000 

RA.7000 

46.50UU 

COWOUCTION 

s 

BETWEEN 

R 

AMD 

10 

23.900D 

82.2000 

" 4. 4000 

CONDUCTION 

2 

BETWEEN 

lA 

AND 

17 

5.0000 

8A.7000 

43*0000 

CONDUCTION 

2 

BETWEEN 

17 

AND 

18 

5.0000 

8A.7000 

35.0000 

CONDUCTION 

2 

BETWEEN 

10 

AND 

19 

A. 0000 

70.0000 

30 .0000 

CONDUCTION 

2 

aETuECN 

10 

AND 

11 

4»ooao 

99*0000 

2G«DCUU 

CONDUCTION 

2 

BETWEEN 

11 

AND 

12 

5.0000 

113.0000 

2A.OOOO 

CONO'jCTION • 

3 

BETWEEN 

Ifi 

AND 

17' 

■ i .0000 ■ 

0.1337 

lA.OOOO 

FORCED CONVECTION 

22 

3ETJEEH 

1 

AND 

20 

55.0000 

36.5000 


FORCED CONVECTION 

22 

HE T’JETfN 

2 

AND 

20 

69*0000 

38.6000 


FORCED CONVECTION 

22 

BETWEEN 

3 

AND 

■ 20 

51.0000 

65.8000 


FORCED CONVECTION 

22 

BETWEEN 

13 

AND 

20 

51.0000 

86.3000 


FORCED CONVECTION 

22 

BETWEEN 

15 

AND 

20 

51.0000 

77.0000 


FORCEO CONVECTION 

22 

BETWEEN 

10 

A NO 

20 

85.0000 

27.0000 


FORCED CONVECTION 

22 

BETWEEN 

11 

AMD 

20 

99.00QO 

28*6000 


FORCED CONVECTION 

22 

BEIWEEN 

12 

and 

20 

109.0000 

35.0000 


FORCED CONVECTION 

22 

between 

lA 

AND 

20 

82.2000 

AO.OOOO 



"o 0“ 

.“ini 


-Q-- 

c<r*>« 

'2. 


O 
O . 

-w-r^ 

«o 

> o 

(— irn . 
-•p*4'TSi‘ 
.03! 
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*•* CYBEAN/MAS A_» TECHNOLOGY DIVISION SKF _ INDUSTRIESJNC. *_ C Y_8,E N /_N ASA 
DESCRIPTION OF THE GEOMETRY AND INDICATION OF THE TYPES AND PATHS OF HEAT TRANSFER 

all lengths are in HILLIHETERS^ a "NEGATIVE SIGN OF THE INDEX'HEANS NO ROTATIONAL" ‘SYHHETRT 

TYPE OF heat TR. INDEX NODE NODE 1ST LENGTH' 2ND LENGTH " '3RD LENGTH" 


FORCED 

CONVECTION 

22 

BETWEEN 

17 

AND 

20 

82.2000 

24.0000 

FORCED 

CONVECTION 

22 

BETWEEN 

■18"' 

ANO 

20 

75.0000 

19.7000 

FORCED 

CONVECTION 

22 

BETWEEN 

19 

AND 

20 

34.5000 

G9.0Q00 

FORCED 

CONVECTION 

23 

BETWEEN 

10 

AND 

30 

87.2000 

23.0000 

FORCED 

CONVECTION 

23 

BETWEEN 

11 

ANO 

30 

101.0000 

41,7000 

FORCED 

CONVECTION 

23 

BETWEEN 

12 

ANO 

30 

115.0000 

35.0000 

FORCED 

CONVECTION 

23 

BETWEEN 

14 

ANO 

30 

87.2000 

40.0000 

FORCED 

CONVECTION 

23 

BETWEEN 

17 

ANO 

30 

87.2000 

50.0000 

FORCED 

CONVECTION 

23 

BETWEEN 

18 

AND 

30 

78.0000 

29.2000 

FORCED 

CONVECTION 

23 

BETWEEN 

19 

ANO 

30 

34.5000 

G9.0000 

FORCED 

CONVECTION 

21 

BETWEEN 

4 

ANO 

22 

GG.OOOO 

18.3000 

FORCED 

CONVECTION 

21 

BETWEEN 

5 

AND 

22 

G4.0000 

15.5000 

FORCED 

CONVECTION 

21 

BETWEEN 

G 

ANO 

22 

G4.0D00 

15.5000 

FORCED 

CONVECTION 

21 

BETWEEN 

7 

and 

22 

13.0000 

241.0000 

forced 

CONVECTION 

21 

BETWEEN 

8 

AND 

22 

78.0000 

27.9000 

FORCED 

CONVECTION 

21 

BETWEEN 

9 

ANO 

22 

72.0000 

33.9000 

FORCED 

CONVECTION 

21 

between 

IG 

AND 

2G 

71.0000 

154.0000 

FLUID 

FLOW 

A2 

FROM 

29 

TO 

21 

(INDEX 41) 


fluid 

flow 

A2 

from 

29 

TO 

24 

(INDEX 41) 


FLUID 

FLOW 

A1 

FROM 

21 

TO 

22 

(INDEX 41) 


FLUID 

FLOW 

Hi 

FROM 

22 

TO 

23 

(INDEX 41) 


fluid 

flow 

A1 

FROM 

23 

TO 

28 

(INDEX 41) 


FLUID 

FLOW 

A1 

FROM 

■24 

TO 

25 

(INDEX 41) 


FLUID 

FLOW 

A1 

FROM 

25 

TO 

2G 

(INDEX 41) 


FLUID 

FLOW 

41 

FROM 

2G 

TO 

27 

(INDEX 41) 
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-C'Y B E A- N N ’A S ' A ’ TECHNOLOGY 0 1 VI SION ' "SKF fNOUSTRIES "INC; * r ■Y'B'r'A’TTTTr'f'S'Ti 

DESCRIPTION OF THE GEOMETRY AND INDICATION OF THE TYPES AN_D_PATHS OF HEAT TRANSFER 

ALL LENGTHS ARE IN MILLIMETERS* A NEGATIVE SIGN OF THE INDEX MEANS NO ROTATIONAL SYMMETRY 


TYPE OF HEAT TR. 

INDEX 


NODE 


NODE- 

1ST LENGTH 

2ND LENGTH 

3RD length 

fluid flow 

A1 

FROM ■ 

27 

TO 

28 

(INDEX 41) 



FLUID FLOW 

02 

FROM 

2B 

TO 

25 

(INDEX 42) 



BEARING CONDUCTION 

51 

SETUCEN 

4 

AND 

7 

1.0000 

1.0000 

l.OOOD 

BEARING CONDUCTION 

51 

BETWEEN 

7 

AND 

8 

1.0000 

2.0000 ' 

■1.0000" 

BEARING CONOUCTION 

51 

BETWEEN 

5 

AND 

7 

1.0000 

5.0000 

1.0000 

BEARING conduction 

51 

BETWEEN 

G 

AND 

7 

l.OOOQ 

6.0000 

1.0000 


oo 

® ■ 
O Si 
O > 
S? sr-. 

so 

r" pi 
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*• * C Y 8 e A'N / N“A S A’ « “ TECHNOCOGY "D IVl SI ON ' ' SKF INDUSTRIi:S‘lMC““ C'VB't "A-W /“N-A'S'TT 


1 

TEMPERATURE MAP 

TEMPERATURES ARE IN DEGREES C_ELS_IUS_. __ T^E_ F IRST 28 TEMPERATURES ARE CALCULATEO; THE OTHERS AR E KNOWN ~ 
STEADY STATE TEMPERATURE C ALCULA TION._ INITl AL TEMPERATURES 

calculated temperatures 


NODE 

TEMPERATURE 

node 

temperature 

NODE 

TEMPERATURE 

NODE 

TEMPERATURE 

NODE 

TEMPERATURE 

1 

93.000 

■ 2 

“ 120.000 

3 

93.000' 

~ 9 

' • 125.000 

• — s — 

125.000'“ 

6 

125.000 

7 

125.000 

8 

125.000 

9 

125.000 

10 

93.000 

11 

. 93.000 

12 

93.000 

13 

93.000 

19 

93.000 

IS 

93.000 

16 

125.000 

17 

93.000 

18 

93.000 ■ 

■ 19 

93.000 ■ 

20 

93.000 

21 

93.000 

22 

93.000 

23 

93.000 

29 

93.000 

25 

93.000 

26 

93.000 

27 

93.000 

28 

93.000 





KNOWN boundary TEMPERATURES 








NODE 

TEMPERATURE 

NODE 

TEMPERATURE 

NODE 

TEMPERATURE 

NODE 

TEMPERATURE ' 

NODE 

TEMPERATURE 


29 93.000 30 29.000 
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ROLLER profile PLOT DESCRIPTION 

ROLLER SLICE RADII ARE PLOTTED AS A 
"FUNCTION OF DISTANCE ALONG' THE ’ROLLER' JC^ATTI S' 

PLOT IS FOR ROLLER GEOMETRY GROUPl 

LENGTHS SKOUN ARE IN MM*S 

SL'ICE RAOII ARE SHOWN BY AN' »S» 

PLOT SHOWN IS ALONG THE ROLLER-OUTER RING EFFECTIVELENGTH 


SCALE MULTIPLIER FOR X-AXIS = O.IOOOCOOE 02 
SCALE MULTIPLIER FOR THIS PROFILE = 0 .999999A E-02 


'C 'O' 

,y »»n* 

o s 

g 5S... 

33 rr 


lO *a_ 

cr> 
> s> 
r- sm 


— 
■< <m 
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^ . „ ■ 

629. S2S 629.810 630.093 630.360 630.665 630.950 631. 23S 631.520 631.805 632.090 632.375 

-0.587 *S ... - - - "4--'- = o-.587 

-O.'iSe + S . ♦ -0.956 

-0.326 + S ♦ -0.326 

-0.196 + S ♦ -0.196 

-0.065 + S ♦ -0.065 

0.065 4 S + 0.065 

0.196'+ '0.196' 

0.326 + S ♦ 0.326 

0.956 ♦ S ♦ 0.956 

0.587 +S ■ + “ 0.587 

4 -4 4 4 — 4 4 4 --- 

SLICE R9DII FOR THIS PROFILE 

629.525' '629Z810' "'630;095 ^63'0';380 630:665 ‘ “eSOigSO 63X7235 ^6317S'20 '6317805 6'327il9o' — ^32;375 ' 
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ROLLER profile PLOT DESCRIPTION 

ROLLER SLICE RADII ARE PLOTTED AS A 
FOfiCTlON OF' distance ALONS "THET ROLLER "X-A'X'tS 

PLOT IS FOR roller GEOMETRY GROUPl 

- 'LENGTHS "SHOUM ARE“IN hH'S'” 

SLICE RADII ARE SHOWN BVAN ‘S*’ 

"PLOT SHOWN IS ALONG THE ROLLER-INNER RING EFFECTIVELENGTH' 

SCALE MULTIPLIliR FOR X-AXIS = O.IOOOOOOE 02 
scale MULTIPLIER FOR THIS PROFILE = 0.999999AE-D2" 
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s ir nrc " gr n rn- TT R 'T m s"t r o n ce 

629.525 6?9.«10 630. D95 630. 3B0 630.665 630.950 631.235 631.520 • 

+ — .i. + — 

-0.5R7 +S 

-O.'iSO + s 

-0.326 + 

-0.196 ♦ 

-0.065 ♦ 

0.065 + 

0.196 + 

0.326 + 

0.956 + S 

0.587 +S 

SLICE RADII FOR THIS PROFILE 

629.525 -629.810 6-30i095 ' 630.380 630.665 ' 6301950 ''631. 235’ ‘ 6311520 “ 


631.805 632.090 632.375 

+ + - + 

' ♦ ' -01587 

+ -0.956 

S ♦ -0.326 

's ♦ ""-o.ige 
S ♦ -0.065 

S ♦ 0.065 

S ♦ ‘ 0.196 

S + 0.326 

+ 0.956 

■ “ 0.587 

+ 1 + 


6311805 6'32709(T ‘ff3'2.375 
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RACEWAY PROFILE PLOT DESCRIPTION 


RACEWAY SLICE RAOII ARE PLOTTED AS A 
FUNCTION OF DISTANCE ALONG THE RING X-AXOS 


PLOT IS FOR ROLLER GEOMETRY GROUP 1 


LENGTHS ARE SHOUN'IM MH*S 


SLICE RADII ARC SHOWN BY AN •S»‘ 


PLOT IS SHOWN along THE ROLLER-OUTER RING EFFECTIVE LENGTH 


scale MULTIPLIER FOR X-AXIS = O.IOOOOOOE 02 
scale MULTIPLIER FOR'THIS PROFILE = 0 ,I 000 B 00 E ‘01 


Sample Output - Page 12 
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•srr c r -k Wf r v » k ' th i 's~r R u k ti n ' 



.5 35 

78.820 

79.105 

79,390 

79.675 

79.960 

80.295 

80.530 

80.815 

81.100 

81.385 


-0.652 

+S 










♦ 

-0.652 

-0.522 

♦s 










4 - 

-0«S22 

-0.291 

♦S 










♦ 

-0.391 

** 0 • ^ & i 

+s 










4 

-0.261 

-O.ISO 

+s 










+ 

-0.130 

-0.000 

♦s 










4 

-0*000 

0.130 

♦s 










4 

'■ 0.130 

0 .261 

♦s 










4 

D«261 

0.391 

♦s 










4 

0.391 

0.522 

+s 










4 

0.522 

0.652 

+s 










4 

0.65'2 





■'SL'rCE 'RAOI t 

FOR THIS 

PROFILE ■■ 







7B.53S 

78.820 

79.105 

79.390 

79.675 

79.960 

80.295 

80.530 

80.815 

81.100 

81.385 
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RACEUAY PROFILE PLOT DESCRIPTION 

RACEWAY SLICE RADII ARE PLOTTED AS A 
FUMCTION of distance' ALONG THE RING X-AXOS" 

PLOT IS FOR ROLLER GEOMETRY GROUP 1 

LENGTHS ARE SHOWN IN HH *S 

SLICE RADII ARE SHOWN BY AN 'S' 

PLOT IS SHOWN ALONG THE ROLLER-INNER RING EFFECTI VE "LENGTH ' 


SCALE MULTIPLIER 
SCALE MULTIPLIER 


FOR X-AXIS = O.IOOOOOOE 02 
FOR THIS PROFILE = O.IOOOOOOE 


01 


Sample Output - Page 14 
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825 66.110 66.395 66.680 66.965 67.250 67.535 67.820 68.105 68.390 68.675 


-0.652 +S 
-0.522 +S 
-0.391 +S 
-0.261 +S 
-0.130 +S 
-0.000 +S 
0.130 +S 
0.261 +S 
0.391 +S 
0.022 +S 
0.652 +S 


+ -O'. 652 

♦ -0.522 

♦ -0.391 

+■■ ‘-o; 26 r 

+ -0.13*0 

♦_ - 0.000 

+ ' 0.130 

♦ 0.261 

+ 0.39i 

+ ■ 0.522 

+ 0.652 


65.825 66.110 


SLICE 'radii FOR'THIS "PROFlCe * — 

66.395 66.680 66.965 67.250 67.535 67.820 68.105 68.390 68.675 


Sample Output Page 15 
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ROLLER TOTAL LENGTH 1A.562 


ROLLER MAXIMUM DIAMETER 12.6A7 


EFFECTIVE L^IGTH (O.Ri> 13.D38" EFFECTIVE LENGTH <I.R.)'' 13.038-' 

FLAT LENGTH 8.AOO CROWN RADIUS 622.300 

ROLLER MASS 0.001176 

END SPHERE RADIUS (LEFT) 381;0DO- •• END SPHERE RADIUS ’(RI 0HT7 ‘SBl'.OO O' 

END PLAY (O.R.) 0.0 EMO PLAY (I.R.) 0.051__ 

DIAMETRAL CLEARENCE 0.127000 

■number of rollers 10 total number OF“ROL'LER-RAC£UAY“SLTtES-“ id — 

- OUTER AND INNER RING DATA 

RACE CURVATURE(O.R.) 0.0 RACE CURVATURE ( I .R. ) 0.0 

PITCH DIAMETER 1AA.363 - ' 

SPECIFIED misalignments D.O (Y-RAOIANS) 0 .0 D 1 ABA (Z-RAOI AHS ) 


FLANGE angle FOR FLANGED INNER RING (DEGREES) 

LEFT SIDE 0.600 RIGHT SIDE -0.600 


outer ring SPEED (RPH) 0.0 INNER RING' SPEED (RPH)- 0 .-lOBSODSAE- 05 ~ 
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c“'r-B"L-A-*i / N « I, A ■» — r n: mu t o ir r i r nrr s r un — sin= n ro i RriviiLy twc: — * — c "T ’ h r s ' w 7" h ' a fl 

CAOr OATA 

ryPE-IUNER RING LAND RIDING 

CAGE POCKET CLEARENCE GEOMETRY GROUP 1 0.221 

CAGE POCKET CLEARENCE GEOMETRVGROUP 2 0.0 

CAGE WEIGHT 0.0 

RAIL LAND DIAMETER 157.952 SINGLE RAIL WIDTH A. 572 

rail land 'DIAMETRAL CLEARENCE 0 .982600 
"CAGE POCKET coefficient OF FRICTION! 0.0700 


MATERIAL PROP ERTIES 


MODULUS OF elasticity 

POISSON’S RATIO 

COEF. OF THERMAL EXPANSI'ON'" 

DENSITY 


“ CAGE ■ 'o.R':; iTR^ RTr; preG'; — 

_0.20£ 06_ 0.20E 0_6 0.20E 0_6 D.20E 06 _O.20E 0^ 

0.30 0.30 0.30 0.30 0.30 

""DTiaE-OA 0.12E-09"' 0'.T2E-0A~ a.l2E-0'T'0'; 12E-0A ' 
0.78E 01 0.78E 01 0.78E 01 0.78E 01 0.78E 01 


FRICTION OATA 

replenishment layer THIv^«<.oa 

OUTER RACE .203E 01 
O.R. FLANGE 


INNER RACE .203E 01 
I.R. FLANGE 0.4> 


NASA LIMITING FRICTION COEFFICIENT 0.070 
NASA lube film thickness FACTOR 50.000 
FRACTION OF LUBE IN BEARING CAVITY ” 0.050 


0i 

ti 
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CT'U'"t"'A "N' / N S' UI'VTSION" "'SlirF' 


— * — rTTrrvTT-Tnrrr 


L-IO FATISUE LIVES(HRS) 

OUTER RING ‘ 0.13E 03 
INNER RING 0.76E 02 
SINGLE ROW BEARING ' 0.52E 02 
BEARING FATIGUE LIFE 0.52E 02 


LUBE life reduction FACTORS 
0»R< UR- 

0.620 . 0.526 


USER INPUT LIFE MULTIPLIERS 
D.R. UR. 

l.OOD' ■ — 1.000 


FILM THICKNESS TO SURFACE ROUGHNESS ’RAT 10 ' 
FOR HEAVYEST loaded ROLLING ELEMENT 


O.R. I.R. 

' 1.7B1 l.AAA 
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ROLLER RACEWAY CONTACT LOADS AT THE INNER RING 


roller 

F-X 

F-Y 

F-Z 

H'X 

H-Y 

H=Z 

1 

O.IAOE 02 

0.356E OR 

0.690E 01 

-0.357E 02 

-0.162E 02 

O.SBSE OR 

2 

-0.616E 00 

0.650E 03“ 

0.A82E 01 

-0.239E 02 ' 

-0.9'26E 01 

■r.757E~03" 

3 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

A 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

5 

0.0 

0.0 

0.0 

0.0 

' 0.0 

0.0 

6 

0.0 

0.0 

0.0 

D.O 

0.0 

0.0 

7 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

8 

0.0 

"0.0 ■ 

0 .0 

■ 0.0 

0.0 

‘0.0 

9 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

IQ 

“0.328E 00 

D.692E 03 

0.3A8E 01 

-0.216E 02 

“0.71UE 01 

0.765E 03 


J-4 

ro 
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UIVISIOR ’SKI= 

“IHUirSTRlES’lNC: 

* 


L A H / W y A — 







ROLLER RACEWAY CONTACT LOADS AT THE 

OUTER 

RING 




ROLLER' 

F-X 


F-Y 


F-Z 


M-X 



■H'Y 

-M-r 

— 

1 

-0.251E 

02 

-O.'i'tOE 

04 

0.6S1E 

oi 

0.4S8E 

02 


0.4G8E 01 

0.320E 

04 

2 

0.236E 

00 

-0.155E 

0 ^ 

0.269E 

01 

0.224E 

02 


•0.1<t8E '(Tl 

=o;783r 

"0 3 

3 

"0*1865 

00 

-0.92EE 

03 

-0.333E 

00 

0.3<llE 

01 


•0.77fE-01 

“0.507E 

01 


-0.187E 

00 

-0.922E 

03 

-0.130E 

01 

•0.269E 

01 


-0.819E-01 

O.llSC 

01 

5 

-0.116E 

00 

-0.922E 

03 

-0.353E 

00 

0.328E 

01 


•0.379E-01 

0.733E 

00 

6 

0*0 


-0*916E 

03 

-0*135E 

01 

'0*3025 

01 


0*0 

0.0 


7 

O.llOE 

00 

-0.929E 

05 

-0.175E 

00 

0.441E 

01 


0.39JE-01 

-0.705E 

00 

■ 8 

0.186E 

00 

-0.917E 

03 

-0.516E 

00 

0.223E 

01 


0.7b3E-01 

'0.784C 

or 

9 

0.1 87E 

00 

-0.92'JE 

03 

-0.D77E 

00 

0.187E 

01 


0.631E-01 

-0.117E 

01 

10 

-0.106E' 

01 

-0.161E 

0^ 

0«202E 

01 

0.163E: 

02 


-0*3b6C 00 

-0*775E 

03 


H 

tvj 


o o 

■n 2 

tj s 
o s 

O 2s> 

^ r“ 

O 2 
c: > 

I” 


s 
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*"),~»"'-rTVT-n"'T}-r "n'x yw"* — r T r Tff rTr m n^ Tr T Ti T?m " ' ~ sir F"“ '"'r ywv i y i V ' it.V '’ iN[. ' r ' " « » r w* ' s ~ « '~ 

INNER RINC APPLIED FORCEStHOH£NTS«DISPLACEHENTS 


AXIAL LOAD = 0.0 


RADIAL 

LOAO(T> = 

0.4AA820E OR 

RADIAL 

HOMENT(Y) = 

0.0 

RADIAL 

HISALIGNHEnTIY) = 

0«0 

RADIAL 

LOAO<Z> = 

Q.O 

RADIAL 

HOMENT(Z) = 

0.0 

RADIAL 

HIJALIGNHENTIZ) = 

0 . a 33333 p -01 


calculated inner ring reactive FORCEStHOHENTSiOISPLACEHENTS 


axial 

LOAD =-. 

3R3E 02 





AXIAL 

TRANSLATION =0’;’0' 



radial 

LOAOtY) 

■=-.R6RE 

04 

RADIAL 

HOHENT(Y) 

=-.105E 02 

RADIAL 

TRANSLATIONIYJ =0.D02_E-01 

radial 

R0TAT10N(Y> 30. fl 

RADIAL 

LOAD(Z) 

=0.750E 

01 

RADIAL 

MOHENT(Z) 

=-.406E 04 

RADIAL 

TRANSLATION(Z) =».323E-03 

radial 

ROTATION(Z) =0i369E-01 


Sample Ou^ut - J’j-ge _21 



■n**T 


TmusTRrcs thc 


T' "ra T A W TT A 5 ' A 


■T^T-K'r-'A-N 7' N — r t r C Tt WgLUGT U ' i V I 51PN Si^ r 


ROLLER END-FLANOE CONTACT LOADS AT THE INNER RIN0< LEFT SIDE 


ROLtCR 

t 

X 

F-Y 

F-2 

M-X 

"«-r - 

- H-Z“ 

1 

O.IIIE 02 

-0.372E 00 

0.720E-01 

0.49SE 00 

0.2//E 02 

0*661C 02 

2 

0.397E 02 • 

0.299E 00 

0.196E 01 

"-0.957C or 

"0.20HE 02 

(r;i9ir*or 

3 

0 • 1D6E 0 0 

0.0 

0.0 

-0.3A3E-01 

0.219E 00 

0«R3K 00 


0.187E 00 

0.0 

0.0 

-0.296E-01 

0.22bE 00 

0*676C 00 

5 

0 •1X6E 00 

0.0 

0.0 

-0.161E-01 

D.1U3E 00 

0*336e'00‘ 

6 

0 • 0 

0.0 

0.0 

0.0 

0 .0 

0 «0 

7 

OtO 

0.0 

0.0 

0.0 

0.0 

0*0 

“3 

0.0 

0.0 

0.0 

D.O 

0.0" 

0.0 

9 

D.O 

0.0 

0.0 

0.0 

0.0 

0.0 

10 

o.iioe 01 

0.257E-01 

O.SAOE-01 

-0.257E DO 

0.2yiE-01 

o.sisf: 01 


TO 

cr\ 


■t3 ^ 

O S 

o B 

50 r- 

^ to 


O Q 

o E 
o 

c E 

•a 
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TP-.,.-ir-r*«-T"T TT — rrc H n uo)i. 7-i m fTsrT0 H -"-S CT nmusTirrtrTJjr^ — * — t, ' " r ' a i a ni"T ' Tj -! T 5 "H ' 


ROULER END-FLAfJGC CONTACT LOADS AT THE INNER RINGt RIGHT SIDE 


ROLLER 

• F-X 

• - F-Y ■ " 

■ F-Z ■ 

H-X 

• • h-y ■ 


1 

O.D 

•♦•• •k Ynip n n 

0.0 

0.0 

0 .0 

0.0 

0.0 

2 

02 

Oi'ilOE'OO"' ■■ 

"~0.19SE 01 

■~-0;609E 01 

' "-OVeiYE 71 — 

-D'ilT^r’OT' 

3 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

5 

0.0 

0.0 

0.0 

0.0 

0.0 

■ 0.0 

6 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

7 

-O.llOE 00 

0.0 

0.0 

-0.205E-01 

-0.9J5E-01 

-0.562E 00 

. g ~ 

-0.186C 00" 

0.0 

"0.0 ■ 

" -0.5<l0E-0r 

■-0.219EOO' 

-0'.629E‘“0 0 

9 

-0.187E 00 

0.0 

0.0 

-0.295E-01 

-C.222E 00 

-0.679E 00 

10 

“0.599E 00 

0.0 

0.293E-01 

-0.917E-01 

-O.OlfiE 00 

-0.185E Oi 


to 

‘ 
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?■ Y W t' A 

TTTnrXTTrT — ILUHNULOUY 
SLIDING SPEEDS AT THE 

ROLLER END-FLANGE’ 

NUUSTRits n«r* • 
CONTACT““' * 

e a n / 

ROLLER 

NO 

FLG 1 

FLG 2 

FLG 3 

FLG 4 

1 

0.0 

0.0 

1290.224 

c.o 

2 

0.0 

0.0 

790.277 

1681.253 

3 

0.0 

0.0 ■ 

’1086.589' 

0.0 

A 

0.0 

0.0 

1508.870 

0.0 

S 

0.0 

0.0 

1783.655 

0.0 

6 ‘ 

O.Q 

0«D 

1917.799"" 

0.0 *■ 

7 

0.0 

0*0 

0.0 

697.070 

8 

0.0 

0.0 

0.0 

1109.983 

9 

0.0 

0.0 

0.0 

1485. 888'" 

10 

0.0 

0.0 

796.975 

1681.320 


LUBRICASIT DATA 

I-* LUBE TYPE’-' "HIL-L-2S6B9* 




TEMP 

dens 

VISCOSITY 

VISCOSITY 

PRESS. VIS. 


O.R. 


147. 

0.912 

■ "(CSTKl 
.24SE 

01 

’ {CPOIS)""' 
.223E 01 

" ’ coef; 

.196E 00 


I.R. 


153. 

0.908 

- .228E 

01 

.207E 

01 

.192E 00 


SULK 


133. 

0.923 

.294E 

01 

.271E 

01 

.211E 00 

o ^ 

J9> 

FLG. 

1 

147. 

0.912 

.24SE 

01 

.2235 

01 

.19RC 00 

-»r-— 
O "0 

> o 

FLG. 

2 

147. 

0.912 

.2'(5E 

01 

.223E 

01 

•198E 00 

FLG* 

3 

145. 

0.914 

.251E 

01 

.229E 

01 

.2005 00 

C ^ 

FLG. 


1^2. 

0.916 

.260C 

01 

.238E 

01 

*202E 00 



Sample Output Page 24 
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n *■ * ■■■ r ‘ y"B~ ‘nr iTT inrr - 1 . ‘ * — rtC Hwqt T ) 'G ^ *i y m5mr'*5 K ~p 


Tramy r mu ’ i Nr. — * — rr-B' Txii ; ‘ W"yT"A 


ROLLING element SPE EDS tNON-CONTACT FORCES tOEFLECTI QMS 


ROLLER 

OUTER 

RING 

DEFLECTION 

ROLLER 

CENTRIFUGAL 

FORCE 

ROLLER 

ORBITAL 

SPEED 

roller 

ROTATIONAL 

SPEED 

NORMAL 

FORCE 

"CAGE POCKET FORCE'S" 
Y-FRICTION 
FORCE 

Y-FRICTION 

moment 

1 

0 .0 

0.927E 

03 

0.9103E 

OA 

-0.1134E 

06 

-0.114E 

02 

-O.BOIE 

00 

0.0 

2 

0.0 

0.91AE 

03 

0.9039E 

04 

•«0.U2f.E 

06 

-0.650E 

01 

D.455E 

00 

0.0 

3 

0.0 

0.925E 

03 

0.9096F 

04 

-0.H3IE 

06 

0.35IE 

01 

D.246E 

00’ 

0.0 

0 

0.0 

0.916E 

03 

0.9052E 

04 

-0.U24E 

06 

0 .403E 

01 

0.282E 

00 

0.0 

5 

0.0 

0.926E 

03 

0.9098E 

04 

-0.1132E 

06 

0.288E 

01 

U.202E 

00 

0.0 

' 6 

0.0 

0.915E 

03 ' 

■•o.boaoe 

04 

■■-0.U23E 

06' 

0.394E 

01’ 

U.276E 

00 

0.0 

7 

0.0 

0.926E 

03 

0.9102E 

04 

-0.1132E 

06 

0.26RE 

01 

O.IRHE 

00 

0»0 

8 

0 .0 

Q.915E 

03 

0.90A7E 

04 

-0.112SE 

0& 

0.288E 

01 

0.202E 

00 

0.0 

9 ' 

0 .0 

0.926E 

03 ■ 

0.9101E 

04 

-0.1131E 

0 6 

0.549E 

01 

0.244E 

00 

0.0' 

10 

0 .0 

0.917E 

03 

0.9057E 

04 

-0.1127E 

06 

-0.329E 

01 

U.231E 

00 

0.0 


’calculated CAGE SPEED SOTS;’ 
CAGE DRIVING TORQUE 2A2.1 
CAGE LAND NORMAL FORCE 0.0 
EPICYCLIC CAGE SPEE0.912E 00 ' 
EPICYCLIC ROLLER SPEED. H3E 06 
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Tr-r 


c y n e k K ' r w v ' 5 ’ R ' t "'~r erHwn. T )GV ' mvis i BN — szr 


INDUS FR Its INC. — * — rTTrTTrTr7TrT~5~)r 





HERTZ 

CONTACT 

STRESS AT 

ROLLER-RACE 

AND ROLLER 

END-FLANGE CONTACT 



ROLLER 

I. RACE 


O.RACE 


.0 

1 

FLO 2 

FLG 3 

— — 

FLG 4 

1 

0.201E 

04 

0.187E 

04 

0.0 


1.0 

0.370E 

02 

0.0 

2 

0.93AE 

03 

0.U9E 

04 

0.0 


1.0 

O.S70E 

02 

0.568E 02 

3 

•0.0 • 


'0.800E' 

"03 

■D.0~ 


:.0' 

"o.saoE 

•or 

0.0 

4 

0.0 


0.79 7E 

03 

0.0 


).0 

0.922E 

01 

0.0 

6 

0.0 


0.797E 

03 

0.0 


1.0 

0.7M3E 

01 

0.0 

•6 

0.0 


0.794E 

03 

0.0' 


D.O ■ 

0.13OC 

or 

o.o"' 

7 

0.0 


0.800E 

03 

0.0 


0.0 

0.0 


Q.769E 01 

e 

0.0 


0.797E 

03 

0.0 


0.0 

0.0 


0.920E 01 

-g- 

0.0 


0.800E 

03' 

0.0 ■ 


0.0 ■ 

“0.0 ■ 


■'0'.921E-0T 

10 

0.952E 

03 

0.121E 

04 

0.0 


0.0 

0.1 S8E 

02 

0.137E 02 


w 

o 
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g§ 
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m 



T*-* — r'T-B-T"A-Tr'rTr'Ti’''5-'''Ji' i ' "iLi;HMOL'i5tT"uiVTSior-*SKr nro omnrg'T vc ; — s — ryHT ' rN ’ V 





lubricant 

film THICKNESS 



ROLLER 

I.RACE 

O.RACe 

FLG 1 

FLG 2 

FLG 3 

FLG 

1 

0.367C 00 

0.4S2C 00 

0.0 

0.0 

0.0 

0.0 

R 

0.575E 00 

0.600E 00 

0.0 

0.0 

0.0 

0.0 

3 

'0.0 

'0.662E~00'' 

0.0“ 

'0.0 

“DVO" 

■o;o' 

4 

0.0 

O.OOOE 00 

0.0 

0.0 

0.0 

0.0 

5 

0.0 

O.063E 00 

0.0 

0.0 

0.0 

0.0 

(, 

0.0 

0.060E 00 

0.0 

0.0 

0.0 

0.0 

7 

0.0 

0.6&3C 00 

0.0 

0.0 

0 .0 

0.0 

0 

0.0 

O.06OE DO 

0.0 

0.0 

0.0 

0.0 

■0 

0.0 

■ 0.662E 00 ■ 

0.0 

'0.0 

0.0 

0.0* 

10 

0.573E 00 

O.590E 00 

0.0 

0.0 

0.0 

0.0 


04 


o o 

-n ?o 


■0 
o 
o . 

33 HT 


Q- 

§ 


lO 

c 


■T0 


5=" 9 , 

r* m 
» 
K 
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TT-I — IT Y '«■ f A"N'7 N"H "S '■«■■< — rtCHNULDGT UIVISIOH SKF 


T T croynmy in r ;”"x — c T’g' T 'n- Ti r ' w r ra ' 


ROLLER 


ROLLER 

SKEW AND TILT 

(RADIANS) 


ROLLER SKEW 

ROLLER TILT 

ROLLER SKEW 

nOLLtn TILT" 

RELATIVE 

relative 

absolute 

AasOLUTf 

-.6773E-0Z 

-.2379E-02 

-.6773E-02 


J.777AE-03 

-.S107E-D3 

0.1532E-02 

0,5S60E-03 

-.395AE-03 

-,i(A57E-03 

0.9878E-03 

_ 

-.AB58E-03 

0.4527E-D3 

0.977AE-03 


-.aS60E-D3 

0.11758-02 

0.5938E-03 



0.H55E-02 

— 


0.23fl2E-03 

0.1179E-02 

=.5167E-03 

— — 

O.A079E-03 

0.9'(36E-03 

-.9753E-03 


0.39A1E-03 

-.A515E-03 

-.9891E-03 

— 

'.301RE-03 

■-.SlOAE-03 

■-.1157E-02 

"D;B'662E:-oy"‘“ 
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BEARING HEAT GENERATION RATES 

SUM OF ROLLING ELT/O.R. CONTACT HEAT GEN. RATES 11A3. 

SUM OF ROLLING ELT/I.R. CONTACT HEAT GEN. RATES 339. 

SUM OF ROLLING ELT DRAG HEAT GEN. RATES 1R77. 

SUN OF rolling ELT/CAGE PKT. HEAT GEN. RATES 23A, 

CAGE RAIL/RING LANO HEAT GEN. RATE 277. 

SUM OF ROLLING ELT/FLC. 1 HEAT GEN. RATES 0. 

SUM OF ROLLING ELT/FLG. 2 HEAT GEN. RATES 0. 

SUM OF ROLLING ELT/FLG. 3 HEAT GEN. RATES 57. 

SUH OF ROLLING ELT/FLG. A HEAT GEN. RATES 89. 
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*** CfBEAN/NASA* TECHNOLOGY OIVtSION SKF INDUSTRIES INC. » CYBEAN/NASA 
TEMPERATURE HAP 

TEMPERATURES ARE IN DEGREES CELSIUS. THE FIRST. 28 TEMPERATURES ARE CALCULATED* THE OTHEW ARE KNOWN 
STEADY STATE TEMPERATURE CALCULATION* FINAL RESULT AFTER A ITERATIONS 

calculated TEMPERATURES 


NODE 

temperature 

node 

TEMPERATURE 

NODE 

TEMPERATURE 

NODE 

TEMPERATURE 

MODE 

temperature 

1 

268.720 

2 

283.358 ’ 

3 

161.780 

■■ 4 ■ 

152.296' * 

5"' 

' 144.489 ' 

6 

1A0.370 

7 

1A6.038 

8 

148.280 

9 

154.854 

10 

146.871 

11 

1A7.376 

12 

150.206 

13 

195.432 

14 

161.887 

15 

238.166 

18 

110. 8A7 

IT 

18A.817 

18 

168.344 

19 

162,429 

20 

188.430 

21 

93.000 

22 

132. All 

23 

171.821 

24 

93.000 

25 

93.000 

28 

10A.769 

27 

116.537 

28 

144.179 






KNOWN BOUNDARY TEMPERATURES 

MODE TEMPERATURE NODE “ TEMPER ATURE" NODE ’TEMPERATURE ’NODE ’"TEMPERATURE "■ NODE"TEHPERATURI 
29 93.000 30 2A.000 
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